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1 Introduction 

A TMD factorization [1-3] generalizes the usual concept of parton density by allowing PDFs 
to depend on intrinsic transverse momenta in addition to the usual longitudinal momentum 
fraction variable. These transverse-momentum dependent parton distributions (also called 
unintegrated parton distributions) are widely used in the analysis of semi-inclusive processes 
like semi-inclusive deep inelastic scattering (SIDIS) or dijet production in hadron-hadron 
collisions (for a review, see Ref. [3]). However, the analysis of TMD evolution in these cases 
is mostly restricted to the evolution of quark TMDs, whereas at high collider energies the 
majority of produced particles will be small-x gluons. In this case one has to understand 
the transition between non-linear dynamics at small x and presumably linear evolution of 
gluon TMDs at intermediate x. 

The study of the transition between the low-x and moderate-x TMDs is complexihed by 
the fact that there are two non-equivalent dehnitions of gluon TMDs in small-x and “medium 
x” communities. In the small-x literature the Weizsacker-Williams (WW) unintegrated 
gluon distribution [5] is dehned in terms of the matrix element 

Y,ix{p\D^UU\zx)\X){X\Dm\Qi_)\p) ( 1 . 1 ) 

between target states (typically protons). Here tr is a color trace in the fundamental 
representation, '^x denotes the sum over full set of hadronic states and is a Wilson-line 
operator - inhnite gauge link ordered along the light-like line 

U{zx) = [oon + zx,-oon + zxi [x,y] = fdu {x-yrA^{ux+{i-u)y) ^^^ 2 ) 
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and D^U{z±) = d'^U{z±) — iA^{oon + z±)U{z±) + iA^{—oon + z±)U{z±). In the spirit of 
rapidity factorization, Bjorken x enters this expression as a rapidity cutoff for Wilson-line 
operators. Roughly speaking, each gluon emitted by Wilson line has rapidity restricted 
from above by Inx^. 

One can rewrite the above matrix element (up to some trivial factor) in the form 

asV{xB,z±) = - n)^^ / +un)\X){X\T'^^{0)\p) (1.3) 

where 


J^^{z±+un) = [oon + z±,un + z±\°''^n^gF'l^^{un + z±) 

F^{zx+un) = n^gF'^^{un + z±)[un + z±,oon + zxY^^" (1.4) 

and define the “WW unintegrated gluon distribution” 

'D{xB,k±) = Jd^zx F^’'’'''>^'D{xb,z±) (1.5) 

(Here {k,z)± denotes the scalar product in 2-dim transverse Euclidean space.) It should 
be noted that since Wilson lines in Eq. (1.1) are renorm-invariant asF(xB,k±) does not 
depend on the renormalization scale /r. 

On the other hand, at moderate xb the unintegrated gluon distribution is defined as 

[ 6 ] 

'D{xB,k±,g) = j(fz± F^^’'''^^V{xB,z±,r]), (1.6) 

asVixB,zx,g) = (z^ + rrn)|X)(X|.F<(0)|p) 

where \p) is an unpolarized target with momentum p (typically proton). There are more 
involved definitions with Eq. (1.6) multiplied by some Wilson-line factors [3, 4] following 
from CSS factorization [7] but we will discuss the “primordial” TMD (1.6). The Bjorken 
X is now introduced explicitly in the definition of gluon TMD. However, because light-like 
Wilson lines exhibit rapidity divergencies, we need a separate cutoff g (not necessarily equal 
to Inxfi) for the rapidity of the gluons emitted by Wilson lines. In addition, the matrix 
elements (1.6) may have double-logarithmic contributions of the type (a^r/Inx^)” while the 
WW distribution (1.3) has only single-log terms (a^ Inx^)” described by the BK evolution 
[8, 9]. 

In the present paper we study the connection between rapidity evolution of WW TMD 
(1.3) at low Xb and (1.6) at moderate xb ~ 1. We will assume k‘j_ > few GeV^ so that 
we can use perturbative QCD (but otherwise k± is arbitrary and can be of order of s as 
in the DGLAP evolution). In this kinematic region we will vary Bjorken xb and look how 
non-linear evolution at small x transforms into linear evolution at moderate xb- It should 
be noted that at least at moderate xb gluon TMDs mix with the quark ones. In this paper 
we disregard this mixing leaving the calculation of full matrix for future publications. (For 
the study of quark TMDs in the low-x region see recent preprint [10].) 
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In addition, we will present the evolution equation for the fragmentation function 

V^{PF,k^,r]) = (1.7) 

asV\f3F,ZF,rj) = jdu ^(0|-Ff + un)\p + X){p + X\F<mQ) 

X 

where p is the momentum of the registered hadron. It turns out to be free of non-linear 
terms, at least in the leading log approximation. 

It should be emphasized that we consider gluon TMDs with Wilson links going to -|-cx) 
in the longitudinal direction relevant for SIDIS [11]. Note that in the leading order SIDIS is 
determined solely by quark TMDs but beyond that the gluon TMDs should be taken into 
account, especially for the description of various processes at future EIC collider (see e.g. 
the report [12]). 

It is worth noting that another gluon TMD with links going to —oo arises in the study 
of processes with exclusive particle production (like Drell-Yan or Higgs production), see for 
example the discussion in Ref. [13]. We plan to study it in future publications. 

The paper is organized as follows. In Sec. 2 we remind the general logic of rapidity 
factorization and rapidity evolution. In Sec. 3 we derive the evolution equation of gluon 
TMD in the light-cone (DGLAP) limit. In Sec. 4 we calculate the Lipatov vertex of the 
gluon production by the J-f operator and the so-called virtual corrections. The final TMD 
evolution equation for all xb and transverse momenta is presented in Sec. 5 and in Sec. 6 
we discuss the DGLAP, BK and Sudakov limits of our equation. In Sec. 7 we demonstrate 
that the linearized evolution equation for unintegrated gluon distribution interpolates be¬ 
tween BFKL and DGLAP equations. In Sec. 8 we present the evolution equations for 
fragmentation TMD and Sec. 9 contains conclusions and outlook. The necessary formulas 
for propagators near the light cone and in the shock-wave background can be found in 
Appendices. 


2 Rapidity factorization and evolution 


In the spirit of high-energy OPE, the rapidity of the gluons is restricted from above by 
the “rapidity divide” tj separating the impact factor and the matrix element so the proper 
definition of Ux is ^ 


= 

Al{x) 


Pexp 


ig / du p^A^Jupi + x_l) 


d^k 


— e{e^-\a\)e-^^-^A,{k) 


( 2 . 1 ) 


where the Sudakov variable a is defined as usual, k = api+ldp 2 +k±. We define the light-like 

2 

vectors pi and p 2 such that pi = n and p 2 = p—^n, where p is the momentum of the target 

^Alternatively, with the leading-log accuracy one can take the Wilson line slightly off the light cone, see 
Ref. [3]. To pave the way for future NLO calculation we prefer the “rigid cutoff” Eq. (2.1) which was used 
for the NLO calculations in the low-a; case [14]. 
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particle of mass m. We use metric = (1, —Ij —I, —1) so p-g = {apfiq+aqf3p)^ — {p^q)^. 
For the coordinates we use the notations x, = and x* = XpP 2 related to the light-cone 
coordinates by x* = y^x+ and x, = y^X-. It is convenient to define Fourier transform 
of the operator 

^ IJ dz* ([oo, z.]r9F:::iz.,z±)r m 

where the index p denotes the rapidity cutoff (2.1) for all gluon fields in this operator. 
Here we introduced the “Bjorken /3s” to have similar formulas for the DIS and annihilation 
matrix elements (/3s = xs in DIS and /3s = /3s = ~ for fragmentation functions). Also, 
hereafter we use the notation [oo, z^,]z = [oo*pi -|- z±, -|- z±] where [x, y] stands for the 

straight-line gauge link connecting points x and y as defined in Eq. (1.2). Our convention 
is that the Latin Lorentz indices always correspond to transverse coordinates while Greek 
Lorentz indices are four-dimensional. 

Similarly, we define 

FriPB^zx) = -Jdz, e-^^‘^^*g{F::iiz,,zx)[z.,oo]ry (2.3) 

in the complex-conjugate part of the amplitude. 

In this notations the unintegrated gluon TMD (1.6) can be represented as 

(p|/r(/3s, ^±).F-^(/3s, Ox)|p + ^P2) = j;(p|-F“^(/3s, z^)|X)(A|^-^(/3s, 0^)|p + 

X 

= - ‘iTT‘^S{0j3Bg^V{PB,z±,p) (2.4) 

Hereafter we use a short-hand notation 

{p\di...dmOi...On\p') = Y.(p\T{Oi...d„^}\x){x\T{Oi...OnW) (2.5) 

X 

where tilde on the operators in the l.h.s. of this formula stands as a reminder that they 
should be inverse time ordered as indicated by inverse-time ordering T in the r.h.s. of the 
above equation. 

As discussed e.g. in Ref. [15], such martix element can be represented by a double 
functional integral 

{di...dmOi...On) 

= jDAD^Di) e-*^QCD(d,^) j (2.6) 

with the boundary condition A(x, t = oo) = A(x, t = oo) (and similarly for quark fields) 
reflecting the sum over all intermediate states X. Due to this condition, the matrix element 
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(2.4) can be made gauge-invariant by connecting the endpoints of Wilson lines at infinity 
with the gauge link ^ 

{p\Fi{l3B,XA_)[x± + oopi,y±+oopi]F‘^\(3'B,yA_)\p') (2.7) 

This gauge link is important if we use the light-like gauge = 0 for calculations [16], 

but in all other gauges it can be neglected. We will not write it down explicitly but will 
always assume it in our formulas. 

We will study the rapidity evolution of the operator 


-^"(/3s,xx)^"(/3s,y±) 


( 2 . 8 ) 


Matrix elements of this operator between unpolarized hadrons can be parametrized as [6] 


e*(^’^)-L(p|jr“’?(/3^,2;_L)J'“^(/3B,0_L)|p + |p2) = 2'K‘^5{i)l3Bg^Tlij{l3B,kr,y) 

'k ‘ k^ \ 

Fij{j3B,k±-,p) = - gijV{j3B,k±,y) + {-^ + gij^jn{PB,k±,y) (2.9) 

where m is the mass of the target hadron (typically proton). The reason we study the 
evolution of the operator (2.8) with non-convoluted indices i and j is that, as we shall see 
below, the rapidity evolution mixes functions D and T-L. It should be also noted that our 
final equation for the evolution of the operator (2.8) is applicable for polarized targets as 
well. 

We shall also study the evolution of fragmentation functions defined by “fragmentation 
matrix elements” (1.7) of the operator (2.8). If the polarization of the fragmentation hadron 
is not registered, this matrix element can be parametrized similarly to Eq. (2.9) (cf. Ref. 
[ 6 ]) 



^(0|^r(-/3ir, zx)\p + X){p + ip 2 + X| J-f (-/3^, 0x)|0) 


= 27rH{0PFg‘^ 


gijV^{/3F,k±,g) + 


2kikj 

m? 


+ Sij 


kl 


'H^il3F,k±,p) 


( 2 . 10 ) 


Note that Pf should be greater than 1 in this equation, otherwise the cross section vanishes. 
As to matrix element (2.4), it can be defined with either sign of Pb but the deep inelastic 
scattering corresponds to /3s = > 0. In our calculations we will consider /3s > 0 for 

simplicity and perform the trivial analytic continuation to negative ,0s in the final formula 
(5.2). 

In the spirit of rapidity factorization, in order to find the evolution of the operator (2.8) 
with respect to rapidity cutoff g (see Eq. (2.1)) one should integrate in the matrix element 
(2.4) over gluons and quarks with rapidities g > Y > g' and temporarily “freeze” fields 
with Y < g' to be integrated over later. (For a review, see Refs. [17, 18].) In this case, we 
obtain functional integral of Eq. (2.6) type over fields with g > Y > g' in the “external” 

^Similarly, this gauge link is implied in Eq. (1.1) which is Eq. (2.7) at Pb = 0. 
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fields with Y < r]'. In terms of Sudakov variables we integrate over gluons with a between 
a = and a' = and, in the leading order, only the diagrams with gluon emissions are 
relevant - the quark diagrams will enter as loops at the next-to-leading (NLO) level. 

To make connections with parton model we will have in mind the frame where target’s 
velocity is large and call the small a fields by the name “fast fields” and large a fields by 
“slow” fields. Of course, “fast” vs “slow” depends on frame but we will stick to naming 
fields as they appear in the projectile’s frame. (Note that in Ref. [8] the terminology is 
opposite, as appears in the target’s frame). As discussed in Ref. [8], the interaction of 
“slow” gluons of large a with “fast” fields of small a is described by eikonal gauge factors 
and the integration over slow fields results in Feynman diagrams in the background of fast 
fields which form a thin shock wave due to Lorentz contraction. However, in Ref. [8] (as 
well as in all small-x literature) it was assumed that the characteristic transverse momenta 
of fast and slow fields are of the same order of magnitude. For our present purposes we 
need to relax this condition and consider cases where the transverse momenta of fast and 
slow fields do differ. In this case, we need to rethink the shock-wave approach. 

Let us figure out how the relative longitudinal size of fast and slow fields depends on 
their transverse momenta. The typical longitudinal size of fast fields is u* ~ ^ where l±_ is 
the characteristic scale of transverse momenta of fast fields. The typical distances traveled 
by slow gluons are ~ where k±_ is the characteristic scale of transverse momenta of slow 
fields. Effectively, the targe-a gluons propagate in the external field of the small-a shock 
wave, except the case <C k\ which should be treated separately since the “shock wave” 
is not necessarily thin in this case. Fortunately, when <C k\ one can use the light- 
cone expansion of slow fields and leave at the leading order only the light-ray operators of 
the leading twist. We will use the combination of shock-wave and light-cone expansions 
and write the interpolating formulas which describe the leading-order contributions in both 
cases. 


3 Evolution kernel in the light-cone limit 

As we discussed above, we will obtain the evolution kernel in two separate cases: the “shock 
wave” case when the characteristic transverse momenta of the background gluon (or quark) 
fields l±_ are of the order of typical momentum of emitted gluon A:_l and the “light cone” case 
when <^k\. It is convenient to start with the light-cone situation and consider the one- 
loop evolution of the operator T'“^(/3 b, x_l)T'“*'^(/3b, ?/_l) in the case when the background 
fields are soft so we can use the expansion of propagators in external fields near the light 
cone [19]. 

In the leading order there is only one “quantum” gluon and we get the typical diagrams 
of Fig. 1 type. One sees that the evolution kernel consist of two parts: “real” part with the 
emission of a real gluon and a “virtual” part without such emission. The “real” production 
part of the kernel can be obtained as a square of a Lipatov vertex - the amplitude of the 
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(a) 


(b) 


Figure 1. Typical diagrams for production (a) and virtual (b) contributions to the evolution kernel. 
The dashed lines denote gauge links. 


emission of a real gluon by the Wilson-line operator 

(/5s, 

= - eTri/3B,x^)A^{k)){lim k^A^P{k)^fi/3B,y±))y"^' (3.1) 

Jo-/ la fc2-S>0 ^ fc2_5.o ■> 

(T^X) 

Hereafter we use the space-saving notation = j 2 ^- 

3.1 Lipatov vertex 

As we mentioned, the production (“real”) part of the kernel corresponds to square of Lipatov 
vertex describing the emission of a gluon by the operator The Lipatov vertex is defined 
as 


Lyiik,y^,/3B) = /lim fc2(r{A“(fc) y^)}) 

A : 2—^0 


(3.2) 


(To simplify our notations, we will often omit label rj for the rapidity cutoff (2.1) but it will 
be always assumed when not displayed). 

We will use the background-Feynman gauge. The three corresponding diagrams are 
shown in Fig. 2. 


k 



k 



k 



Figure 2. Lipatov vertex of gluon emission. 


3.1.1 Emission of soft gluon near the light cone 

In accordance with general background-field formalism we separate the gluon field into the 
“classical” background part and “quantum” part 

Af^ -b 
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where the “classical” helds are fast (a < a') and “quantum” helds are slow (a > a'). It should 
be emphasized that our “classical” held does not satisfy the equation = 0; rather, 

= —gTp'yyF'ij;, where ijj ^re the “classical” (i.e. fast) quark helds. In addition, in 
this Section it is assumed that the slow helds are hard and the fast helds are soft so one 
can use the light-cone expansion. We will perform calculations in the background-Feynman 
gauge, where the gluon propagator is ^^7 see Appendix A. 

The hrst-order term in the expansion of the operator [ 00 , y*]y^F^{y^:, y±) in quantum 
helds has the form 


[oo,y*]r^.T(y*,y±) = ~[^,y.];^Ar{y.,y±) 


(3.3) 


roo 2 

- [oo,y,]^^d,Ar{y*,y±) + iJ d-zi [ 00 , y^)^, y^) 


(to save space, we omit the label from classical helds). The corresponding vertex of gluon 
emission is given by 


hm A:2(A“^(A:)([ 00 ,y±)) 


lst\ 


(3.4) 


= lim 


fc2-S>0 


d 


- —[oo,y,]^-(A“^(fe)Ar(y„y^)) - [ 00 ,y,]^-(A“^(fe)5,Ar(y*,y±)) 




dz^ [ 00 , 2:*]j;(A“‘1(A;)A2(z*, y±))[z*, y*]y™F™(y*, y_L) 


y* 


To calculate the r.h.s. we can use formulas (10.47)-(10.48) from Appendix A. As we men¬ 
tioned, we need contributions to production part of the kernel with the collinear twist up to 
two. However, it is easy to see that the light-cone expansion of gluon emission vertex starts 
with the operators of twist one (~ F,i) since the gauge links in the hrst term in Eq. (10.20) 
cancel in Eq. (3.4) and the remaining background-free emission of gluon is proportional to 
<5(/d_B + which vanishes for /3b > 0. Thus, to get the contribution to the production part 
of the kernel of collinear twist up to two it is sufficient to use formula (10.20) for Feynman 
amplitude and formula (10.23) for complex conjugate amplitude with twist-one (one F,i) 
accuracy. In this case the quark terms do not contribute and the gluon terms simplify to 

hm A:2(A“^(fe)A^^(y)) = - y,, fe), 

Ci^i/(oo, y*, y_L;/c) = y^y[oo, y=i,]j; 

+ g dz* i - ;^k^{z - y)*gy„[oo,z*]yF,j{z*,y±)[z*,y^]y 

Jy, ^ 

4 • • \ 

+ -^{dlP2u - dlp2y)[oo,z^]yF,j{z^,yp)[z^,y^]y j 


- 8 - 


(3.5) 





With the help of this formula Eq. (3.4) reduces to 
hm y±))''') 


(3.6) 


+ ^^P 2 ii){[oo,y*]yF,i{y^,y^)[y^,oo\y 


as 


- -j dz^[oo,z^]yF,i{z^,yj_)[z^,oo\y\ 

as ./ y^ / 


25 

+ — 
as 


- df^ki) j dz^ {[oo,z^]yF,j{z^,y±)[z^,oo]y)'""''^ 




Note that k^ x (r.h.s. of Eq. (3.6))^ = 0 as required by gauge invariance. Integrating the 
r.h.s. of Eq. (3.6) over y* we obtain 


Lt(k,yi,M = i^lim i^(.47(t)(J-?(/3B.!/i))“> = 


(3,7) 


ajdss / A:^ \ d d , 

- api^)h - 5^ki + 


J^r{PB + ^,y±) 

as 


I I ajdBsgyik’- 2 akik’- 

- + hi + ah. + q+«&-/''■ J 

At this point it is convenient to switch to the light-like gauge = 0. Since ky x 

(r.h.s. of Eq. (3.7))^ = 0 it is sufficient to replace ap^ in the r.h.s. of Eq. (3.7) by 

ap^ — k^ = — k^ — ^P 2 - One obtains 

(3.8) 

\kySl d^ykj + S\k^ - Qyjk^ klgyjk^ + 2 k^kik^ ^ 


L A)2 


+ —,y±) + 0 {p 2 y) 

as 


a/IfiS +(a/3ss + fc^)2 

We do not write down the terms ~ p 2 y since they do not contribute to the production 
kernel (~ square of the expression in the r.h.s. of Eq. (3.8)). 

For the complex conjugate amplitude one obtains from Eq. (10.49) 

fe2 

-i lim k‘^{Al{x)Al{k)) = A:), 

k^—^O ^ ^ 

/ Ai 

Oyy{x^,oo,x±;k) = gyAx*,oo]x + g dz^ [x^, z^j^ci—^{z - x)^gyuF,j{z^,x±)[z^,oo]xk^ 

. / /r> ' G~S 


4 . . ^ 

- ■^i^iP 2 u - dip 2 y)F,j{z^,X±)[z^, 00 ], 


(3.9) 


where Oyi, is obtained from the Eq. (10.23) with twist-two accuracy (as we mentioned, 
quark operators start from twist two and therefore do not contribute to the production 
kernel). 

Repeating steps which lead us to Eq. (3.8) we obtain 

= -i lim = 2gF^^'^^^ (3.10) 

^■2 — 


ky5^ S^kj + A:^(jf - gyjk^ k\gyik^ + 2kji[kik^ - 

■ k\ af5BS + k\ {oifiBS + k^Y 


+ + 0 {p2y) 
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The product of Lipatov vertices (3.8) and (3.10) integrated according to Eq. (3.1) gives the 
production part of the evolution kernel in the light-cone limit. To get the full kernel, we 
need to add the virtual contribution coming from diagrams of Fig. lb type. 


3.2 Virtual part of the kernel 

To get the virtual part coming from diagrams of Fig. lb type we need to expand the 
operator T up to the second order in quantum held 




(3.11) 


4^2 /-oo /-z* 




'y* 


'y* 


8 

+/ dz^{[oo,z^]y{A'^{z^,y±)[z^,y^]y)''"'A'l^'^{y^,y±)) 
jy* 

^ Jy* 

As we mentioned above, we are interested in operators up to (collinear) twist one. Looking 
at the explicit expressions for propagators in Appendix A it is easy to see that the only 
contribution of twist one comes from (A 2 ( 2 *, y_L)A3(y*, y^)) propagator, which is given by 
Eq. (10.45) with 


Q,i{z^,y*-,pi_) = - —[ dz'^ [z^,z'^]F,i{z'^)[z'^,y^] 

Jy. 


(3.12) 


We obtain 


([oo,y*]--F.7(y„y^))2-d = — - (y±|^|y±)[oo,y.]--F.7(y„y^) 


da 


1 


J V* 


0 “ L 


P± 


dz. (y±|e-4(-*-?/*)|y^)[oo,z,]-™T.7(z„y^)J (3.13) 

where we used Schwinger’s notations 

{x±\f{p±)\y±) = Jd'^px edP^^~y'>^f{p), ( x _ l | p ±) = 

For the operator T{j3B,yi.) the Eq. (3.13) gives 


(3.14) 


(.^r(/?B,y±))'"" = -g^Nj —{y^\ 


da , 


af^BS 


0 a p]_{al3BS + p\) 


For the complex conjugate amplitude 

1 mn\ 2nd 
/ 

r*oo 


^2 I'oo 


\y±)J7il3B,y±) (3.15) 


(3.16) 


= - ^/ dz^ I dz'^ F^{x^,x±){[x^,z^]x{A^{z^,x±)[z^,z'^]xA^{z'^,x±))[z'^,oo]x)"" 

^ ■I'y-t J Zt 

8 f°° ~ 

- * 5 ^ J dz*{A'^^{x^,x±){[x^,z^]xA'^{z^:,x±)[z^,oo]x)"^"') 

9j r°° 

H- g dz^{diA^'^[[x^,z^]xA’^{z^,x±)[z^,oo]x)"^"') 

^ J Xt, 
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Again, the only contribution of twist one comes from a:^)) given by Eq. 

(10.46) with 

Qii{x^,z^]PA_) =- — f dz'^ {[x^,z'^]F,i{z'^)[z'^,z^]f^ (3.17) 

J IC* 

(see Eq. (10.23)) so the virtual correction in the complex conjugate amplitude is propor¬ 
tional to 

W(/3b, 2^±)>""'* = -/jv./”—mi— 

Jo a p^j_{al3BS + pj_) 

The total virtual correction is 


Wb , x± )Fj {/3 b, y ±))virt 

= - 2g‘^N,FnPB,x^)F^{l3B,y±) T— [ 2 , 2 , (3-19) 

Jo ot J Pl(apss+Pl) 


Note that with our rapidity cutoff in a (Eq. (2.1)) the contribution (3.19) coming from the 
diagram in Fig. lb is UV finite. Indeed, regularizing the IR divergence with a small gluon 
mass m? we obtain 



a/3BS 

(p^ -|- m?){al3BS + p\+ rrF) 





al3BS + m? 
m? 


(3.20) 


which is finite without any UV regulator (the IR divergence is canceled with the corre¬ 
sponding term in the real correction, see Eq. (3.24) below). This feature - simultaneous 
regularization of UV and rapidity divergence - is a consequence of our specific choice of 
cutoff in rapidity. For a different rapidity cutoff we may have the UV divergence in the 
remaining integrals which has to be regulated with suitable UV cutoff (for example, see 
Refs. [20, 21]). Let us illustrate this using the example of the Fig. lb diagram calculated 
above. Technically, we calculated the loop integral in this diagram 

J ad[3d/3 ^ P± _j_ _ iFj{a(3s — -|- ie){aj3's — p]_ — wJ + ie) 

(3.21) 


by taking residues in the integrals over Sudakov variables f3 and f3' and cutting the obtained 
integral over a from above by the cutoff (2.1). Instead, let us take the residue over a: 




(f2p L 


m- 


= */3_b 


+ P1, 


a/3d:(3' 


e{/3)e{-/3') - 9{-/3)e{/3') 

{f3> + l3B-ie){/3-ie)i/3'-f3) 
dpd/3' ^e{/3)6{-(3')-0{-/3)e{/3') 


(3.22) 




+ p]_J /3'+ I3b - iel {/3 - ie){/3'- P) 


+ 


OiP') 


{f3-ie){/3'-P + ie)\ 


[ d^p± f d/3di3' _ 9{/3) ^ r d^pi_ d/3 

J m?+p\J (3'+/3 b - ie{/3 - ie){/3'- 13+ ie) ^JrrF+p^Jo l3{/3 +/3 b) 


which is integral (3.20) with the replacement of variable /3 



as ‘ 
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A conventional way of rewriting this integral in the framework of collinear factorization 
approach is 


r (t^p± r°° di3 
Jm‘^+p]_Jo /3{/3 + Pb) 


f dz 

J m? +p\Jq 1- z 


(3.23) 


where z = is a fraction of momentum {/3b + /d)p 2 of “incoming gluon” (described by /Fi 

in our formalism) carried by the emitted “particle” with fraction /3bP2, see the discussion of 
the DGLAP kernel in the next Section. Now, if we cut the rapidity of the emitted gluon by 
cutoff in fraction of momentum 2 , we would still have the UV divergent expression which 
must be regulated by a suitable UV cutoff. 


3.3 Evolution kernel in the light-cone limit 

Summing the product of Lipatov vertices (3.8) and (3.10) (integrated according to Eq. (3.1)) 
and the virtual correction (3.19) we obtain the one-loop evolution kernel in the light-cone 
approximation 

'-a^k2. + ’^,X2.)/Ft{/3B + ^,y±) 

V G/S O/S 

+ 5^kik^ + 5\kjk^ - 5\k 


2ivt dta 
= 2g^Nj —c 
.L' a 




25^5] 


L 


-k 



(3.24) 

as 


— d^kjF — g^^kikj 

- Qijk^k'’ 


a/dBS + k"^^ {a/dBS + k^Y 

2g,jk'^k^ + 5^^kjk^ + 5\kik^ - 6pik^ - 6lkjk^ k\g,jk^k^ 


+ ki 


{a/dBS + k\)^ 
a/^BS 


3Fi {13 b, X2 _){/3B,y± 


{a(3BS + 

I Ina^ 


k‘j_{a(3BS + k'j_ 

where rapidities of gluons in the operators in the r.h.s. are restricted from above by Incik 
Let us write down now the evolution equation for gluon TMDs defined by the matrix 
element (2.9). If we define /3b as a fraction of the momentum p of the original hadron we 
have /3b < 3. Moreover, in the production part of the amplitude we have a kinematical 
restriction that the sum of /3 b and the fraction carried by emitted gluon ^ should be less 
than one. This leads to the upper cutoff in the k± integral k'^ < a(l — //b)s and we get 
the equation 

^ :{p\Tn/3B,x^)T-{/3B,y±)\p) 

d‘^k± {ed’^'^~y'^^{p\F^{/3B + ^,x±)/Ft{PB + ^,y±)\p) 


din a 


(3.25) 


TT 

5^5^- 

^ 3 


25^6] 


+ 


k\5^j\ + 6^kiF + dlkjk’^ - 5%k^ - SfkjF - g’^^hkj - gijFF 


+ k 


k\ (t/3bs + kj_ {(J/Sbs + k\Y 

2gijFF + 6^kjF + 5]hk^ - - 5[kjk^ k^gijk^F 


{a/3BS -k )3 


0{l - (3 b- —) 
as' 


a^BS 


kj_{a//BS + k^ 


{a/3BS + k\Y. 
-{p\3'i{l3B,X2_)F‘'{/3B,yi. 
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(there is obviously no restriction on k± in the virtual diagram). 

If the target hadron is unpolarized one can use the parametrization (2.9) 


{p\^i (Pb, Z± ) J""(/3b, 0_l) Ip + 

= 27 r^' 3 ( 0 / 3 Bff^ - gipWB^z^^ri) 
= 27r^(5(C)/3BC/^[-fi'yfP(/3B,2±,77) 


■^{2didj + gijd\)n{^B,z^,ri) 
-^{2ziZj + gijzj_)'H"{/3B,z±,r]) 


(3.26) 


where p = In cr, (/3b, - 2 ±,p) = fd'^k± e k±, g) and 

V."{13b, z±,g) = (gfi)^'H(/3B, ?/)• Rewriting Eq. (3.25) in terms of variable (3 = ^ one 

obtains 


dg 


gijasV{l3B,z±,g) H- K{2ziZj + gijz]_)aa'H''{^ b, z±,g) 


m 


(3.27) 


a 


,Nr 


TT 


d/3 ^gijJo{\z±\\/cr(3s)as'D{/3B + I3,z±,g) 

3/3 2/32 


/3b + /3 2 3/3 2/32 /33 

PPb 13b Pb{Pb + /3) /3b(/3b + /3)2 ^ /3b(/3b + /3)3J 

/3 


+ J2(|2;±|\/o^)( 2^ +5ij)Q!sI?(/3B + /3,2;±,??)-^^^^ 


+ ^Jo(|2±| yo)^)(2zi2:j +5jj2;i)as?7"(/3B + /3,2;±,??) 

+ ^^zj_J2{\z±\-\/af3s)as'H''{(3B + f3,z±,g) 


13 


13 


2,02 


[3 b + /3 2 

- /3b/3 13b /3b{/3b + I3)- 

/33 


/3b 


d/3 


TT Jo /3(/3 b + /3) I' 


I3b{I3b + (3) /3b(/3b + /3)2 /3b(/3b +/3)^- 


gijas'D{l3B,z±,g) H - 2(2zjZj + 2:_l, ?/) 


m 


where we used the formula 


1 




- J2{kz)k\ 


i\ziZj+gi^ 


(3.28) 


The evolution equation (3.27) can be rewritten as a system of evolution equations for T> 


































and H" functions {z' = 


^as'D{l3B,z±,r]) 


(3.29) 


agNc dz 


TT 


^ dz' 
Hb 


r{-/o(k±iyfTS/3B^^) [(y^)+ + ^ " 2 <^sV{^,z^,r,) 


+ ;^(1 - z')z' z\j 2 {\z±\\j asjdB ^ J ^asH” , z^.-q)^, 


■^asn"{PB,z±,r]) 


OgNc dz 


TT 


+ 


'/3b 
1 — z' 


4zi z‘ 


T{jo[\z±\^ash^^-^) [(y^)+ - l] asn''{^,z^,rj) 

7^J2{\z±\^asf5B^-^)asV{^,z±,r))^ 


where dzf{z)g{z)+ = dzf{z)g{z) — dzf{l)g{z). The above equation is our final 
result for the rapidity evolution of gluon TMDs in the near-light-cone case. 

It is instructive to check that the evolution equation (3.29) agrees with the (one-loop) 
DGLAP kernel. If we take the light-cone limit x_\_ = (44> z± = 0) we get 


A 

dg 


as'D{/3B,0±,g) 



dl_ 


( 


1 

1-z' 


) 



2 + z'(l 





(3.30) 


One immediately recognizes the expression in the square brackets as gluon-gluon DGLAP 
kernel (the term j5(5(l — z') is absent since we consider the gluon light-ray operator mul¬ 
tiplied by an extra as). It should be mentioned, however, that Eq. (3.30) is not a proper 
DGLAP equation since the latter is formulated for the gluon parton density on the light 
cone dehned by 


dg{xBMg^) = - du (3.31) 

where the light-ray gluon operator F^(un)[un,0]F°‘’‘{0) is regularized with counterterms at 
normalization point /i^ (recall that on the light ray T-product of operators coincide with 
the usual product). 

Gomparing Eqs. (2.4) and (3.31) we see that dg{xB) = F{xb,z±_ = 0) modulo dif¬ 
ferent cutoffs: by counterterms for dg{xB,g^) and by “brute force” rapidity cutoff Y < g 
in D{xb,z±_ = 0,?7 = In a). However, with the leading-log accuracy subtracting the coun¬ 
terterms is equivalent to imposing a cutoff in transverse momenta of the emitted gluons 
k\ < g?. If we would calculate the leading-order renorm-group equation for the light-ray 
operator -T“(/3s, x_l)J-'?(/3b) 3/±) we would cut the integral over from above by g? and 
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leave the integration over rapidity (a) unrestricted. Thus, we would obtain 




(3.32) 




IT 


10 


a J tiPi 

as 


I fffi Bn B) 

rl 

2/3s 

h 

0% 1 

|C71^1 PB PJ 


(0B+00 

' (0B + 00 

{0b + 0Y\ 


X (p|TT(/3b + /3,xx)-F™(/3b + /3,^±)Ip)^ 
„2 


Pb 


PiPs + fi) 




= ^nJ 

TT Jo 


da 


dP 


0s 


Bn B) 

rl 

2/3s 

, 0l 

0B 1 

PB p) 


i0B+0r 

i0B+00 

{0B+00i 


X {p\^{Pb + /3,x^)T^\/3b + /3,x^)\p)>^ - 


0B 


0i0B + 0) 

which should be compared to Eq. (3.24) with x±_ = y_L 


{p\^?{0B,Xx)J^^\l5B.Xx)\pr'] 


{p\K{0B,Xx)r^{0B,Xx)W'^^ 

= 


TT 


T-/ a 




-1 


X {p\^{0B + 0, Xx)J^^00B + 0, X±)|p)'""' 


Wb 0l 0B 1 

Wb + pt- ^ W+W ^ W+W- 


2 

In the leading log approximation 0 ~ 0b = xb so one can replace the cutoff ^ in Eq. (3.32) 

by the cutoff = a and hence dg(xB,ln/x^) = D{xb,z± = 0, Inds) with the leading-log 
accuracy. The equation (3.32) can be rewritten as an evolution equation 


P^^{p\^i0B,Xx)r^{0B,Xx)\pr 

Ots{p) 


(3.33) 


TT 


iVe rd0 {0(1 -0B-0) 


1 


2/3b 


-0 {0B 

X {p\F0{0b+ 0,xx)J^\0B+ 0,xx)\pr - 


0? 

0B 


+ 


0 


0 


0{0B+0) 


{0B+0? {0B+00\ 

{p\J0{0B,Xx)r'-{0B,Xx)\pr} 


which can be transformed to the standard DGLAP form [22] 


2 d 


-^as{p)dg{xB0'a 0^) 


^ JdB 


( 


1 

\-z' 


)+ 



2 + z'{l-z') 


as{0)dg{^M0^) 


(3.34) 


There is a subtle point in comparison of our rapidity evolution of light-ray operators to the 
conventional p? evolution described by renorm-group equations: the self-energy diagrams 
are not regulated by our rapidity cutoff so the (5-function terms in our version of the DGLAP 
equations are absent. ^ Indeed, in our analysis we do not change the UV treatment of the 
theory, we just define the Wilson-line (or light-ray) operators by the requirement that gluons 

®For Eq (3.34) the absence of these terms is accidental, due to an extra ois in the definition (2.4). 
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emitted by those operators have rapidity cutoff (2.1). The UV divergences in self-energy 
and other internal loop diagrams appearing in higher-order calculations are absorbed in 
the usual Z-factors. So, in a way, we will have two evolution equations for our operators: 
the trivial evolution described by anomalous dimensions of corresponding gluon (or 
quark) fields and the rapidity evolution. Combined together, the two should describe the 
evolution of DIS structure functions. Presumably, the argument of coupling constant in 
LO equation (3.30) (which is by default) will be replaced by a^BS in accordance with 
common lore that this argument is determined by characteristic transverse momenta. ^ We 
plan to return to this point in the future NLO analysis. 

4 Evolution kernel in the general case 

In this section we will find the leading-order rapidity evolution of gluon operator (2.8) 

with the rapidity cutoff Y < r] = In a for all emitted gluons. As we mentioned in the 
Introduction, in order to find the evolution kernel we need to integrate over slow gluons with 
a > a > a' and temporarily freeze fast fields with a < cr' to be integrated over later. To this 
end we need the one-loop diagrams in the fast background fields with arbitrary transverse 
momenta. In the previous section we have found the evolution kernel in background fields 
with transverse momenta l± <C p± where p± is a characteristic momentum of our quantum 
slow fields. In this section at first we will find the Lipatov vertex and virtual correction 
for the case l± ~ p± and then write down general formulas which are correct in the whole 
region of the transverse momentum. 

The key observation is that for transverse momenta of quantum and background field 
of the same order we can use the shock-wave approximation developed for small-x physics. 
To find the evolution kernel we consider the operator (2.8) in the background of external 
field x±) (the absence of x, in the argument corresponds to a = 0). Moreover, we 

assume that the background field A, (x*, x±) has a narrow support and vanishes outside 
the [—CT*, CT*] interval. This is obviously not the most general form of the external field, but 
it turns out that after obtaining the kernel of the evolution equation it is easy to restore 
the result for any background field by insertion of gauge links at ±oopi, see the discussion 
after Eq. (5.4). 

Since the typical d’s of the external field are Bfast ~ ^ ^ ^ the support of the shock 
wave u* is of order of ~ This is to be compared to the typical scale of slow fields 
g ^ ~ ^ ^ cj* SO we see that the fast background field can be approximated by a narrow 

shock wave. In the “pure” low-x case = 0 one can assume that the support of this shock 
wave is infinitely narrow. As we shall see below, in our case of arbitrary we need to 

^Note that while in the usual renorm-group DGLAP the argument of coupling constant is a part of LO 
equation, with our cutoff this argument can be determined only at the NLO level, same as in the case of 
NLO BK equation at low x [14]. This is not surprising since we use the rapidity cutoff borrowed from the 
NLO BK analysis. 
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look inside the shock wave so we will separate all integrals over longitudinal distances z* in 
parts “inside the shock wave” |z*| < fi* and “outside the shock wave" | 2 ;*| > u*, calculate 
them separately and check that the sum of “inside” and “outside” contributions does not 
depend on cr* with our accuracy. 


4.1 Production part of the evolution kernel 

In the leading order there is only one extra gluon and we get the typical diagrams of Fig. 
3 type. The production part of the kernel can be obtained as a square of a Lipatov vertex 



Figure 3. Typical diagrams for production (a) and virtual (b) contributions to the evolution kernel. 
The shaded area denotes shock wave of background fast fields. 

- the amplitude of the emission of a real gluon by the operator T"? (see Eq. (3.1)) 

where the Lipatov vertices of gluon emission are defined as 

i“-(fc,y±,/3B) = i lim k\A';,{k):F^{/3B,y±)) 

LY^{k,x±,l3B) = -i lim (^s,x_L)i“(fe)) (4.2) 

(cf. Eqs. (3.2) and (3.10)). Hereafter {O) means the average of operator O in the shock- 
wave background. 

4.2 Lipatov vertex of gluon emission in the shock wave background 

As we discussed above, we calculate the diagrams in the background of a shock wave of 
width ~ ^ where l±_ is the characteristic transverse momentum of the external shock-wave 

field. Note that the factor in the exponent in the definition of J-{(5b) is ~ /3b^ which is 
not necessarily small at various /3b and and therefore we need to take into account the 
diagram in Fig. 4c with emission point inside the shock wave. We will do this in a following 
way: we assume that all of the shock wave is contained within cr* > z* > —cr*, calculate 
diagrams in Fig. 4a-d and check that the dependence on cr* cancels in the final result for 
the sum of these diagrams. 
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Figure 4. Lipatov vertex of gluon emission. 


We start the calculation with the expansion of the gluon helds in J-{I3b, z±) in the hrst 
order in slow “quantum” held: 


•^r(/5s,y±) == 9 I d-y^ e 


JpBV* 


2ig 


- + — / dz^ [oo,z^]yA^{z^,y_L)[z^,y^]^'^F^{y^,y_L) 


'y* 


(4.3) 


where the gauge links and are made of fast “external” helds. The corresponding vertex 
of gluon emission is given by 

hm e{A-\k)FnPB,y±)^^^) = lim k^g [dy, - i/3s[oo, 

- -[oo,y,]--5,^r(2/*,y±) + ^ rdz, ([oo,y±)}> (4.4) 


The diagrams in Fig. 4a, 4b, and 4(c-d) correspond to different regions of integration over 
y* in Eq. (4.3): y* > n*, — cr* > y*, and cr* > y* > — cr*, respectively. 

The trivial calculation of Fig. 4a contribution yields 


i lim /c^(^“^(A:)T7(/3B,y±))Fig.4a 


fc2^.0 


(4.5) 


= ig dy,F^^y*[oo,yX^\ime{A-{k)[-i/3BAT{y.,y±)--diA:^{y.,y±)]) 


= gd 


( 7 * 

an 


9^liCt|dBS + 2akiPiy_ ^i{l3B+*y^)o-,-i{k,y)j 

ot/dBS + k\ 


4.2.1 Diagram in Fig. 4b 

Next step is the calculation of Fig. 4b contribution. Using the vertex of gluon emission 
from the shock wave (11.30) one obtains 

ighm [ "‘dy, 

k^^QJ-oc S 

= 9 [ dy^ [(fz± - iPBOyi{oo,y^, z±; k){z±\F^y*\y±) 

J —oo J 


2i 


- -0^.((X),y„zx;fe)(zx|y,e4?/*|y^)}“-(C/t)™n 
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where O is given by Eqs. (11.31): 


0 °^y{sx),y^,zx_-,k) g^^U'. 

+ g f dz^ ([(X),z*] 2 [ - - iD^) + -^{Sj^p 2 L, - 6 ip 2 iM)]F,j{z^, z±)[z^, -oo]^)"^ 

/_Ol5 G^S 


ab 

z 


+ 


4g 


dz^ ([ cx ), 2 *] 2 { - ip2ii.P2vD^ F,j{z^,z±)[z^,-oo\2 


+ 9 

^2 ^ /*oo 


/ ^* 4 h 

dz'^ [2iagfj,yz'^ - P2^lP2l] F,j{z^,Z2_)[z^,z^]zFj{z'^,Z2_)[z'^,-oo\:,]Y 

-oo ® 

g2 , roc pzt 

+ ^\ dz^ I dz(V'(z*,2;_L)[z*,oo]2t“C/2i^[-oo,4]27M 

^J—oo J—oo 


' —oo J —oo 

roo roo 


dz^ dz'^ip{z^,z±)'y^ ^i'y^[z^,-oo]zt^UlF[oo,z'^]:,ip{z'^,z±)'^ (4.7) 


' —OO J 


where we replaced y* by —oo since we assumed that there is no gauge held outside the 
[—(T*,(T*] interval. 

Let us compare relative size of terms in the r.h.s. of this equation. The leading 
term is ~ t/z ~ 1 and it is clear that all other terms are small. Indeed, the hrst term 
in the second line is ~ f dz.^ z^{2k^ — iD^)F,j ~ ^kWjU ~ ^ 1 since the width 

of the shock wave is ~ p- and a ^ a' (recall that in this Section l± ~ k±). Similarly, the 

hrst term in the fourth line is ~ ^ Jdz^dz'^ z'^F,j{z^)F,^{z'^) ~ ^d^UdjU ~ cr*^ <C 1. 

Next, let us hnd out the relative size of quark terms in Eq. (4.7). The “power counting” 
for external quark helds in comparison to gluon ones is 

^Jdz^ijj ■]^)iil){z^) ~ ^Jdz^D^F,i(z^) ^ k'j_U r-u A:^ and each extra integration inside the 
shock wave brings extra c*. Thus, the two last terms in Eq. (4.7) are ~ S'/ll/<C 1 
After omitting small terms the expression (4.7) reduces to 


4y • . 

0^u{oo,-oo,z±;k) = g^yUz+ — 2 ^^ 1 P 2 u - dlp 2 ^) dz^ [oo,z^]zF,j{z^,z±_)[z^,-oo \2 

J—oo 

Ag f 

H- y^P 2 plP 2 u dz^ [(X), 2 ;*] 2 l - iD^F,j{z^, z±)[z^,-oo]^ 

OL S J— oo ^ 

- y / dz'^ F,j{z^,z±)[z^,z',]zFj{z',,z±)[zi,-oo]z'^ 


= gfzuUz + ^^{6lp2u-6ip2fz)djUz-‘^^^dlUz (4.8) 


as 


where we used the formula 

roo 

d^Uz =g dz^ [oo,2*]j 


(4.9) 


X (^^D^F,j{z^,z±)[z^,-oo]z+ ^ I dzi F,j{z^,z±)[z^,zi]zFj{z'^,z±)[zi,-oo]z'j 


89 


®Note, however, that the quark term ~ ^ 2 ^ 3 P 2 tiP 2 iyJdztD^F,j{z,) ~ P 2 ^,P 2 v is of the same order 


of magnitude as the gluon term ^ 


iP 2 ^P 2 v f dz,dziF,j{z*)F,\zi). 
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Using Eq. (4.8) one obtains for the r.h.s. of Eq. (4.6) 


i lim A:^(A“'^(A:)J’”(/3s,y_L))Fig.4b 


(4.10) 


= g [d^z±e g U* 

J aj3BS+p\ 


Pias 




a/lss + pj_ 


= - ge "^^'"*{k±\[gf,iU - 


as 


diU 


+ + -d^U - ^^dlU)pi] 

Pbs a ^ a‘^s ' -' 


af^BS 


a/3ss +pj_ 


u^\yi.r 


P'2 

where we used the fact that ^<7* <C 1 when all the transverse momenta are of the same 
order. 


4.2.2 Diagrams in Fig. 4 c,d 

Next step is the calculation of Fig. 4 c,d contributions. Using the vertex of gluon emission 
from the shock wave (10.47) and Eqs. (11.6), (11.7) one obtains 

hm kHg rdy, (A“^(fe){ - i/3s[oo, 2 /^) (4.11) 


2. 

-po 

s 


y*]y'^diA'^^{y^,y±)+ ‘^ f dz^ {[oo, z^]yA’^{z^,y±)[z^,y^]y)'^'^F^{y^,y±)}) 

^ J-Ut 


= 9 J dy^ - ‘^{iki +-^)0‘l,'^{oo,y^,y±;k)] 


-iik,y)± 


•5 Jy:, 

= g,g-*(fc,y)x f e''dBy*{^-i{l3^0yi{oo,y^,y±-k)+ -kiOy,{oo,y^,y±-,k)}[y^,oo]y 

J-at L '® 

2 d 4ig f°° 

- -^(C’M.(oo,y*,?/_L;/c)[y*,oo]y) + dz^ 0^,{oo,y^,y±] k)[y^, z^]yF,i{z^,y±) 

4ig f°° 

X [z^joo]^^) + — / dz* 


= +2 — 2* 


^*0^,(oo, z^,y±-, k)[z^,y^]yF,i{y^,y±)[y^,oo]y 


where Oyi, = G^u + Qfiu + Qy.u and G, Q and Q are given by Eqs. (11.6) and (11.7). As 
we mentioned above, the contributions with extra (z — cr)* are small and so are the quark 
terms (except term ~ D^F,j). So, we have Qyi = Qyi = 0 and 


d/ji(oo, ?/*, ?/j_; fc)[y*, oo]jy — g^i 




as^ 


dz^ [oo,z^]yF,i{z^,y^)[z^,oo\y, 


'y* 


/■“ 2 

Oy,{oo,y^,y±-,k)[y^,oo]y = pi^ + g dz4oo,z^]y{ — F,y{z^,y±)[z^,oo]y (4.12) 

Jy» CKS 

- ‘^:^^^D^F,j{z^,y±)[z^,oo]y - f dz'^F,j{z^,y±)[z^, z%F,^{z'^,y±)[z'^,oo]y} 


a^s^ 
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After some algebra the r.h.s. of Eq. (4.11) reduces to 


i lim k^{Al^{k)Fl^{PB,yi_))FigAc+d 


= _ + 2akiPi^ rgi(fe + ^)<7. _ g-i(/3s + ^)cr*1g-i(fc,y)x 

a^BS + k‘j_ ^ 

+ + hP2^.6i - h5l) [T^iUB^VA.) - idjUyUle-^f^^'^*] 

+ - idyUyUy^^^-* - ^[Vil3B,y±) - e-*^«-‘aiC/,C/t]} 

+ ^ fdzM, 

PBS^J 


X [oo, z^]{-^F,y,{z^,y±)[z^, z% - F,j{z^,y±)[z^, z%}F,i{z'^,y±)[z'^,oo]i 


as 


- [dzMjz': 6{z, - 


a^s 




X [z”,z%F,i{zi,y±)[zi,oo]y 


a^BS 


^^^^*{d^UydiUl + 



( 4 . 13 ) 


where 


V(/3i?,y±) = g 



— OO 
OO 


dz^ (^[oo^z^]yD^F,j{z^,y±) 


dz'^ [oo, z%F,j {z'^ , y ±) [z(, z^]yF,^ (z*, y± 


,oo]y ( 4 . 14 ) 
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4.2.3 Lipatov vertex 


The sum of Eqs. (4.5), (4.10), and (4.13) gives the Lipatov vertex of gluon emission in the 
form 


Lfi{k,y±,PB) = ^ - .n ' ' e 


^ ^^^ab dfiiO^PBS + 2akiPl^ 
af^BS + k\ 

+ ^ U^\y±T^ 

'■ as ^ a ^ a^s ' pbs-' apBS + Pj_ 

+ + l3BP2y6l - h5l) [J),(/3s,y±) - idjUyUle-^^^^*] 

- ‘^MPB,y±) + ^^[V{PB,y±) - 

+ ^/!{^f.ifiB,y±) - id^UyUle-^^^^* - ?^[Vi/3B,y±) - 

+ 7^ f dz^dzi e^dBminiz.,z'P 
PbS^J 

X [oo,z^]{—F,y{z^,y±)[z^,z'^]y- 


as 


as 


v2^2 


32ff^P2At 
O^S^jdB . 


dz^dz^dz” 9{z^ — z”)e 


a^s 

'>\f,idBi^in{z'^,z") 


F,j{z^,y±)[z^, zi]y}F,i{zi,y±)[zl,oo\ 


[oo, z^]yF,j{z^,y±)[z^, z%F,\z”, y±) 


X [4,zi]yF,i{zi,y±)[zi,oo]y - 


a^BS 


o-i^B 


^*{d^UydiUl + ^^dlUydiUl)]"'’ (4.15) 


This expression explicitly depends on the cutoff a^. However, we can set cr* = 0 in the 
r.h.s. of Eq. (4.15) (and eliminate few terms as well). To demonstrate this, let us consider 
two cases: /3 r <C — and 13n > — ■ In the first case 

' ^ CT^ ' ^ - CT^ 


Lf,(k,yi,M (fe. yi) + 9{ki I 


-[ gpiU - ^diU + -F(p u + -djiu- ^8it/)K] 

as pbs a ^ ^ 


ajdBS 


a‘^s 


al3BS + p\ 


u^\yz) 


ab 


= 9{k±\gyi{ 


afdBS 


- U ^ 2 apiy[ 


afdBS + p]^ al3BS + p\_ 


- U- 


Pi 


aj3BS + p\_ aj3BS + p\_ 


t/t) 




and all other terms are small since they contain extra factors e^^BZ, _ g-ifeo-, ~ (z — cr)* 
(or — fj* or z” — fJ*) in the integrand. 

In the second case a'(3 b s > p\ so a(3BS 3> p\ and we get 


L%{Ky^,PB) 


— ^-iik,y)A 




+ 


8g^ 

a(^BS^ 


dz^dz( e»femm( 2 *,<) [oo, z^]F,y{z^,y±)[z^, z(]yF,i{z(,y±)[z(,oo]y) (4.17) 


ab 
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where we used the formula 




a^BS 


Pi 


- U 


Pi 


a 2 +2api^{ ^ 2 

'a/3BS+p\ a/3BS+p\ ^ ^a(3BS+pj_ 

= {k± I {2ik^djU - dlU)—^ -+ 2iapif,diU 


a^BS+p^ 


u^)\y±) 

1 


a(3BS + k\ 


aPBS + pj_ 


a/3BS+pj_ 


uMy±) 

(4.18) 


Let us now compare the contributions of various terms in the r.h.s. of Eq. (4.15) to the 
production part of the evolution kernel defined by the square of Lipatov vertices (4.15). It 
is clear that the square of the first term ~ (^ — is proportional to and 

^ p2 ^ 

contributions of all other terms are down by at least one power of . Thus, with our 
accuracy 




o/3ss>p'j^ 


= " 2e 


-i[k,y) 


(^-^)-^“'(/3B,y±) 

as k. 


(4.19) 


We see that in both cases (4.16) and (4.19) one can replace n* by 0. Moreover, with our 
accuracy the Lipatov vertex (4.15) can be reduced to the “direct sum” of Eqs. (4.16) and 
(4.19): 


Lfi{k,y±,l3B) = 2ge 


+ g{k±\gyi{ 


—i{k,y)± _ OtPly 

^ as 

af^BS af3BS 


{^-^mWB,y±)-U,iy±)] 

^ as k, 


ab 


a/3BS+p]_ 


- U 


+ [2ij3BP2ydiU -2idj[Upi^ - 

^ as 


R , r , 2 

apss+p^ apss + pj_ 

2p2p „o n 1 


Pi 


- U- 


Pi 


aj^BS +p\ 


PI 


aj3BS + p]_ 

M rr L. \ah 


(4.20) 

C/t) 


where we introduced the notation Ui = -P'j(O) = i{diU)U'^. It is clear that at /3scr* <C 1 the 
first term in the r.h.s. of this equation disappears and we get the r.h.s. of Eq. (4.16). On 
the other hand, as we saw above, at > 1 all terms in the last two lines in the r.h.s. 

of Eq. (4.20) are small except (fc_L|^^17j — which cancels the second term 

in the first line of Eq. (4.20) so we get the r.h.s. of Eq. (4.19). It is worth noting that at 
(5b = 0 Eq. (4.20) agrees with the Lipatov vertex obtained in Ref. [23]. 

It is instructive to check the Lipatov vertex property k^L'jk{k,y±, (5b) = 0. One 

obtains 


k^ X (r.h.s. of Eq. (4.20))^ 
a(5BS 


= g{k±\ki 


- U- 


a(5BS 


a(5BS + k\ a(5BS+p\ 


-C/t)+fei( 


ki 


-U- 


Pi 


a(5BS + k\ a(5BS + p\ 


+ {a(5Bs\pi, U] + \p\ + a(5BS, U]pi] 


aPBS+Pj_ 


U^-Ui\y^) = 0 


(4.21) 

f/t) 
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4.3 Lipatov vertex for arbitrary transverse momenta 

Let us demonstrate that for arbitrary transverse momenta the Lipatov vertex of gluon 
emission is given by the following “interpolating formula” 




■'fii 


= g{k±\ 

2i 


a^BS + A;^ 


{2ik^djU - dlU) 


H- P2^idiU 2 

as apBS + p\^ 


a^BS 2p2f,. 


- [2idf,U - -^dlU] 


aj3BS+pj_ 
Pi 


[/'*' + 2iapiadiU 


a/^ss +pi 


c/t 


as 


a/Sss +pj_ 


t/t - 


p\ 


+ 


2ge r _ • 2akik^pi^ a(3Bsg^,ik^ 

q;/?bs + ^ * a^BS + k^ a^BS + k‘j_ 


+ |3BP2^J.Sl 




kl 'J 


X [J-,(/3b + (4.22) 

Let us consider at hrst the light-cone limit corresponding to the case when the char¬ 
acteristic transverse momenta of the external “fast” gluon helds are small in comparison 
to the momenta of “slow” gluons which we integrated over. As we discussed above, the 
higher-twist terms DjF,k or F,jF,k exceed our accuracy so we can eliminate terms 
^ dj_U and commute operators djU with ^6®^lting in 


r.h.s. of Eq. (4.22) 


+ 


+ 


P2iiI3b 
ajSBS + k\ 

2ge-^(^^y'> 


If. _ I _ TJ — 

* aj3BS + k\_ 

2akiypi, 


{a/3BS + kj_y 
Ui\y±r 


Ui + 


apifj. 


aj3BS + k\ 


■u,. 


k\ 

al^Bsg^ik^ 


- r Sill , Fi-fi , '.^h'B^y^l^r^- , n 1-7 o( 3 bS ^ji 

■pif, + H [ - + « + a&,, + If + 


X [77(/3B + ^,!,i)-CJ,(9i)]“‘ (4.23) 

It is clear now that the hrst two lines in the r.h.s. cancel the last term in the square brackets 
in the last line so we recover the light-cone result (3.8). 

Next we consider the case when the transverse momenta of fast and slow helds are 
comparable so the Lipatov vertex is given by Eq. (4.20) above. The difference between the 
r.h.s.’s of Eq. (4.22) and Eq. (4.20) is 


2ge r ■ 2akikypif, aPBsg^,ikF ., ^ 

al5BS + k‘^P ^ * al3BS + k\ + A;^ j\ ^ 




(4.24) 


where we used Eq. (4.18). It is easy to see that the expression (4.24) is small in both 
Bb ^ ^ and Bb > ^ cases. Indeed, when Bb ^ ^ the integral representing FABb + 

^,y_L) — Uj{y 2 _) contains an extra factor _ i ^ <C 1 in the 

, ^2 

integrand and in the (3b > — case the Eq. (4.24) is ~ FpfF <C 1 in comparison to the 
leading term in this limit (4.19). 

As in the light-cone case, for calculation of the evolution kernel it is convenient to go 
to the light-like gauge Since A:^ x (r.h.s. of Eq. (4.22))^ = 0 (see Eq. (4.21)) 
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it is sufficient to replace ap^ in the r.h.s. of Eq. (4.22) by ap^ — = — k^ — ^P 2 - One 

obtains 




-like 


a|3Bsg^,i-2kj^ki . ^ 


— -- 2ik-tdiU— -^ 

apBS+Pj_ ^ apBS+Pj_ 


(4.25) 

c/t 


- 2id^U 


Pi 


, 2k'j7 , 

apBS+pj_ kj_ 


+ 2ge 


-i(k,y) 


\kjl5l 
V k^ 


Sfiki + S{k^ - g^ik^ 

a^BS + 


k'igyik^ + 2kj^kik^ 
{uPbs + kPy 


X [J^j(/3B + ^,y±) - C/j(y±)]“^ + 0{p2y) 

As usual, we do not display the term ~ p 2 ^ since it does not contribute to the evolution 
kernel. Using U] = — 2ijPdjU + d‘j_U one can rewrite this vertex as 


Lt{k,y^,PBf^^^-^^^ = i lim A:2(A“q(fe)J-f(/3i,,y^))'i*5ht-iike 


A:2-s>0 


(4.26) 


apBS + p^j_ apBS + k^ 


as 




a^BS+k\ 
where we introduced the notation 


Siki + 6 {- g^ik^ ^ gyik\k^ + 2k^kik^ 


{a^BS + A:^)2 


^f{l3B + —,y±) + 0{p2y) 


as 


- kp k"^ 

J^j {PB + —,y±) = J^j {Pb + —,y±) - (0,y±) 

OtS OtS 


(4.27) 


(recall that Fj{Q,y_i_) = Uj{y±) = idjllylly). 

It should be emphasized that while we constructed the Lipatov vertex (4.22) as a 
formula which interpolates between the light-cone result (3.8) for small transverse momenta 
of background fields and shock-wave result (4.20) for comparable transverse momenta, we 
have just demonstrated that with our leading-log accuracy our final expression (4.22) is 
correct in the whole range of the transverse momenta. 

It is convenient to rewrite the Lipatov vertex (4.26) in a different form without explicit 
subtraction (4.27). Starting from Eq. (4.25) we get 




fll 


(4.28) 


^ as^ I a^BS + kj_ aPBS+pj_ 


Pi 


- 2k^U—^—- 2gyjU ^ „ 

^ aPBS+p]_ aPBS + p]_ 

where the operator J-i{P) is defined as usual 


[/t + M 




9ij 


t/t 

+ 0{p2y 




(4.29) 
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Let us prove that Eq. (4.28) coincides with Eq. (4.26) with our accuracy. First, as 
we discussed above, in the light-cone case {l‘j_ p‘j_) we can drop higher-twist terms and 

commute operators U with pi and ^2 ^ which gives us Eq. (3.8). Second, consider 

the “shock-wave” case ~ When /3 _b W the integral representing 


contains an exponential factor e ~ Ic* 'ptijg factor can be approximated 

k2 

by one, since <C 1 in the shock-wave case (see the discussion above), so we can replace 
+ ’^) by Uj and get 

^ ~ fJ^^\ k,id^u + id^Up,)—^ —^ 

'I a^BS + kl ^^a^BS + pl 

Pi 2ki: 


-U^ - 2kj^idiU 


- 2idkJJ- 


-U^ + 




-u, 


a/3BS + p]_ 


(4.30) 

c/t 


a^BS+p^ 

which gives the first two lines in the r.h.s. of Eq. (4.25). As it was shown above, the last 
two lines in the r.h.s. of Eq. (4.25) are small at 13b <C W sq Eq. (4.28) coincides with 
Eq. (4.25) at Pb 'C W with our accuracy. Finally, in the 13b > ^ case a(3BS and 

therefore the Eq. (4.28) reduces to 


Lf,{k,y^,l3Bf^^^-^^ 


-like 


2k^ 

_^e-^iKy)^gJ:f^|3B + ^^y^) + 0{p2i.) 
K I OtS 


which is the same as Eq. (4.25) in this limit. 

Similar calculation for complex-conjugate amplitude gives 


(4.31) 


= g{xi.\u 


p\9kii + 2pjYPi _ p\gtii + 2pj^pi 


a^BS +p]_ 


k^ 




aPBS+p\ 


|fe±)'“ + + -±,x±) 


as 


- 2ge 
where 


i(k,x)^ + k^6l - g^jki ^ g^jkiy + 2kik^kj^ ^ ^ M 


oPbs + k\ 


{a(3BS + kj_)‘^ 


{Pb + —Fx) + 0(P2m) 
■' as 


~ k'^ -k"^- 

Fj{l3B + —.xa) = Fj[(3B +—,x±) - Fj{Y),XA_) 

CX,S GtS 


(4.32) 


Similarly to Eq. (4.28) we can rewrite the above expression in the form without subtractions 




-like 


(4.33) 


U 


1 


a/^BS +pi 

- 2g^jU 


2 - 2kjYgijU 


1 


af3BS + k\ 


Pi 


a/3BS +pj_ 
k^ 






The production part of the evolution kernel is proportional to the cross section of gluon 
emission given by the product of Eqs. (4.26) and (4.31) integrated according to Eq. (3.1). 
To find the full kernel we should calculate the virtual part. 


26 - 




















4.4 Virtual correction 

To get the virtual correction shown in Fig. 5 we should use the expansion (3.11) of the 
operator T up to the second order in quantum field. From Eq. (3.11) one gets 


(-^”(/?B,2/±))'“" (4.34) 

^OO 

Jv* 


= — [dy, 


2i 


-/ dz^{\oo,z^]y{A^{z^,y±)%,y^]yY'^diA^'^{y^,y±)) 


'y* 


45 


dz^ / dz'^ {[oo, z^]y{A^{z^,y^)[z^, z%A^{z'^,y^)[z'^,y^]y))'^'^F^{y^,y^) 


As in the case of production kernel we will calculate the diagrams in Fig. 5a, 5b, and 5c 
separately and then check that the final result does not depend on the size of the shock 
wave cj* (it is easy to see that the diagram in Fig. 5d vanishes in Feynman gauge). 



(a) (b) (c) (d) 


Figure 5. Virtual part of the evolution kernel. 


4.4.1 Diagram in Fig. 5a 

Let us start with the diagram in Fig. 5a. Using Eq. (3.11) and (11.27) we get 




(4.35) 


= - / dy^ 




poo 

Pb / dz^{{A2{z^,y±)Uyy"'A'P'^{y^,y±)) 
J (J* 


2 ? 1 
- - / dz,{{A^,{z,,y^)UyY^diAr{y,,y^)) 

^ J at 

= f d^ j f dzt{y±\e~"^''*{(^BG,iioo,-oo;p±) 

^ Jo ^ J — OO J (T^ 

+ ^[^••(oo, -oo;p_l) + Q..(oo, -oo;p_L)]7'i}e*^^®+^)*'*f7^'|5_L)^* 

1 .P^\ 

= d:a{y±\^e~^^‘^*{°‘^BsG,i{oo,-oo-,p±) 

Jo PI 

+ 2a[g,,{oo, -oo;p±) + Q,,{oo, -oo]p±)]pi}—~^-~^e~ddB+^y*ut\y_^'^ki 

aPBS + p]_ 
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(as usual we assume that there are no external fields outside [cr*, —u*] interval). Moreover, 
from Eq. (11.25) we see that Q,i{oo, —oo;p±) = — ^diU and from Eqs. (11.25), (11.7) 

and (4.9) that G»»{oo, —oo]p±) + Q,,(oo, —oo;p±) = — ^dj_U so we obtain 

-g [ (4.36) 

^ J — OO 

Jo « P± a/JBS + Pl 

= r — iy±\^[a(JBsd,U-idlUp,] ^ , U^\y±f 

Jo « P± a(JBS+p^j_ 

(recall that -^cr* <C 1 if the transverse momenta in the loop are of order of transverse 
momenta of external fields). 


4.4.2 Diagram in Fig. 5b 

To get the contribution of the diagram in Fig. 5b we need the gluon propagator with one 
point in the shock wave (11.8), which we will rewrite as follows 


/ O ^ r\ 


yj-j! — I „ iva-Li- - (4.37) 

J —OO 

X + Ql%{y*,z^\P±)]\z±) + ( 5 ;±|Q^p(y*, 2 *;p±)e"*^(^"^)‘|y_L)} 

with Q and Q given by Eqs. (11.6) and (11.7) 


as 

Gll{y*,z^]P±) = - 


rz^ 

/ d4([2;*,4]77,i(2;')K,y*])“> Qi,{y*,z*;p±) = Qi,{y*,z^-,p±) = 0, 

Jv^ 


a^s^, 

ig 


dzi / dz'^{[z^, zi]F,j{zi)[zi, z”]F,^{z'^)[z'',y^]) 


ab 


(4.38) 


'y» 


Q,,{y*,z^-PB) = - ^ dz(([2:*,4]D^F,j(4)[z(,y*])“^ Q„{y*, z^-,p_i) = 0 

and therefore from Eq. (4.34) we get 


(4.39) 


25 


2 r—cFt. 


dy^F^^y- Pb dz^{[oo,z^]y{A'^{z^,y±)[z^,y^]y)'^"^A'P'^{y^,y±)) 


- — / dz^{[oo,z^]y{A'^{z^,y±)[z^,y^]y)"''^diA^^{y^,y±)) 

^ J -cr* 

First, let us show that the second term in the r.h.s. of this equation vanishes. From Eq. 
(4.37) we see that 


{A':^{z.,y±)d,At^iy.,y±)) (4.40) 

= -ij ^{y±\Pie~"^‘'^~'''>*[Gt'^{y*,z^-,p±) + Qll{y^,z^;p±)]\y±) = 0 
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because operators in Eq. (4.37) do not contain p and {yA_\pie = 0. 

Now consider the first term in the r.h.s. of Eq. (4.39). From Eq. (4.37) we get 

(^“(z*,y_L)-4-(y*,2/±)) 

and therefore 


^(y±|e [ *dz'^ [z*,z'^]F,i{z'^)[z'^,y^]\y^y'^ 

« Jy, 


(4.41) 


^ J—oo ^ J—a* 


.•i. 


f°°d:a r 

X / — iy±\[e * - 


■fl / 

g as ^ 


/o a 


] 2, |yx) (4.42) 


pi(a/3BS+pi) 


where -7f'«,y±) = fi(([oo, <]yE,i( 2 ;(, y_L)[z(, oo]^)^', see Eq. (1.4). Since |^(t* <C 1 the 
r.h.s. of this equation can be simplified to 




(4.43) 


= -2 — 


s Jq a 


\y±) ( dz[ {z'- a)^Fi\z'^,y±) ~ 0 


a/3ss +pj_ 


because it is 0(—cr*) in comparison to Eq. (4.36). 


4.4.3 Diagram in Fig. 5c 

As in previous Sections, we start from rewriting Eq. (3.11) 




(4.44) 


-O'* 

2 rcr. 


= ^ / dy* e 


l/3sy. 


45 


fZ, 

dz^ / d^ ([oo,z*]j^(A^(z*,y_L)[^;*,4]?/^2(4,y±)) 

Jv* 


fy* Jy* 
roo 


POO 

X [4,5*]y)””'4C(5*,5±)+ /3 b / d2;*([oo,2;*]y(A2(2;*,y±)[2;*,y*]y)’'”'A™"‘(y*,y_L)) 

Jy* 


2i 


-/ dz*([oo,z*]2^(A2(2;*,y_L)[2;*,y*]y)”™5iA™‘i(y*,y_L)) 


'y* 
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Using the propagator (11.8) with point y inside the shock wave (and point z anywhere) ® 
we obtain (hereafter di{C) = —i[pi,C]) 

fcr* _ „ 

iy ^ 91 tJJ-j / (4.45) 

(r — 1 /V,. r n / , 

[oo,z^\ 






+ 


^9 


X G,i{z^,y*-,p±)[y*,oo])^^ + ([oo, z* ] O,, (z*, y*; p_L) [y*, oo]) 

- (9i([oo,2;*])C>,,(z*,y*;p_L)[y*,oo] - [cx), z*]C>„( 2 ;*, y*;p_L)9i([y*, oo])| J|y_L) 

[ (y± I e"*^ ([oo, z*] O,, (z*, ; p_L) [^1, y*](y*) [y*, oo]) Iy_L) I 


where G,^ is given by Eq. (11.6) 

2 g f^* 

G,iix^,y^;p±) =-— / dz*[a:*,z*]F,j(z*)[z*,y*] 

Jy* 


as 


(4.46) 


and O 99 = + Q 99 by Eqs. (11.6), (11.7) 

g r^* 

0..(a;*,y*;p±) = — dz* [x*,z*]{ - iZl'^F,j(z*)[z*,y*] 


a^s 


'y* 


4y 


f dz( F,j(z*)[z*,z(]F,^(z()[z(,y*]} (4.47) 

Jv^ 


Using these expressions, one obtains after some algebra 
^ Hcly, e*^«^‘([oo,y,]"-F.7(y„y^))^'g-5'= 

S .1 _rr 


(4.48) 


= 9^"“ /“-(Sxl 


1 


|y±){ / dw, 


0 « ''''^'p\{a/dBS + p‘i) 

/ 1 2i \ 

X ( - 2zQ;/3Bd^i(w^*,y±) - -d^V{w^,y±) + ^V(u;*,y_L) / d2(l^i(4,y±)) 

^ ^ J Wif 

‘2i%q /2z^ 

— / dz( / drc* ( — [oo,u;*]j;Zl'^F,j(r(;*,y_L)[w^*,oo]y [e 

® y-CT* y^t ^ s 


+ 


_ g-*/3BCT* 


g^2 /-lu* 

H—^ / dr(;([oo,u;*]yE,j(r(;*,y_L)['i«*,w^*]t/E,-^'(t(;()[t(;(,oo]y 


where, as usual, = y([oo,'u;*]yF,j(r(;*, y_L)[t(;*, 00 ]^)^^ and 


kl 


V('w*,y±) = 2iy[cx),t(;*]yiy-^F,j(u;*,y_L)[t(^*,oo]y 

g^2 I'oo 

H-/ du;( [cx),u;(]j/F,j(u;(,y_L)[w'*,w^*]y7^,^(^«*,y±)[^i'*,oo]j^ (4.49) 

^ Jon. 


(cf. Eq. (4.14)). 


^Strictly speaking, one should depict Eq. (4.44) as several diagrams with points z (and z') inside and 
outside the shock wave. 
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4.4.4 The sum of diagrams in Fig. 5 


The total virtual correction coming from Fig. 5 is given by the sum of Eqs. (4.36) and 
(4.48) 




Fig. 5 


(4.50) 


= r — {y^\^[af3BsdiU -idlUpi]-^^ , 

a a(3BS + pj_ 


10 


da 


g 


'i 


U^\y^r 
) 


+ 


+ /r'= 7 ^ 

/ 1 2i \ 

X {-2ial3BFi{w^,y±)--dfV{w^,yA_)+ ^V(u;*,y_L)/ dz'^Fi{z'^,yi_)\ 

S ^ J Wt ^ 

dz[ ! dw^ C^^[oo,w^]yD^F,j{w^,y^)[w^,oo\y\e'^^’^^* — 

J-Ut Jz'^ ^ S 

gg2 1‘Wf 

H- ^ dw'^[oo,w^]yF,j{w^,y±)[w^,w%Fj{w'^)[w'^,oo]y 

^ J z' 

X [e‘/'«<-e-*^«"*])j-,(4,y±)} 

Let us prove that with our accuracy it can be approximated as 


kl 


^5 1 dy*e*^-^*([oo,y,]-™F.7(y„y^))Pig- ^ (4.51) 

= nV”"' r — {y±\^[af3Bsd,U-idlUp,] ^ ^ U^\y±)’^^ 

Jo a al 3 BS + pj_ 

_ ■g2fnki r^(^^| ,A y±m{/3B,y±) - id^UyUlf 

Jo a pj_{af3BS + pj_) ^ 

To this end we compare the size of different terms in the r.h.s. of equations (4.50) and 
(4.51) at (Jb^^* ^ 1 and ^ 1- In the first case (at /3 b(T* <C 1) the only surviving terms 
in the r.h.s.’s of these equations are the hrst terms and they are obviously equal. 

In the second case let us start from Eq. (4.51). Since /3 _b(T* ~ /3b ^ > 1 we have 

P± 

uPbS /§> p\ so 


r.h.s. of Eq. (4.51) = - ig^T^{y^\\\y^)Ff{PB,y±) (4.52) 

Jo « PI 

Let us now compare the size of different terms in the r.h.s. of Eq. (4.50). Since 4 Jdw^V(w^) 
~ d]_l7l7^ the hrst term in the fourth line ~ Jdw^,a(JBF»iiw*,y±) ~ ajJBsdiUyUy is 
much greater than the second term ~ ^fdw^diV{w^,y±) didlUyUl or the third term 
~ ■^fdw^diV{w^,y±)fdw'^F,i{w'^) d‘j_UydiUy. Moreover, it is easy to see that the terms 
in the last three lines in Eq. (4.50) are of the same order as the terms ~ V in the fourth 
line so they are again small in comparison to the term ~ Fi. Thus, we get 


r.h.s. of Eq. (4.50) = pkl^-r0B<r, / _^y^\ ^Q.u_ Q 2 up^p\\y^^kl 

Jo a aPBS 

poo ^ 1 o pa* 

- / —{y±\^\y ±)-/ dw, (u;„y^) (4.53) 

Jo « P± 
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which coincides with the r.h.s of Eq. (4.52). 

Last but not least, let us prove that one can use the formula (4.51) in the light-cone 
limit l‘j_ <C p\ where it coincides with Eq. (3.15). First we notice that the term ~ d^U has 
twist two and so exceeds our twist-one light-cone accuracy. Next, since the commutator 
[p‘^,diU] consists of operators of collinear twist two (or higher), one can rewrite the first 
term in the r.h.s of Eq. (4.51) in the form 


{y±\^aPBsdiU 

PI 


a^BS +pj_ 




iy±\ 


a/^BS 


pj_{al3BS+pi) 


\y±)diUyUl (4.54) 


so it cancels with last term in the r.h.s of Eq. (4.51) and we obtain 




(4.55) 


= - ig f 


2 j^nkl 


a^BS 


’ da 

~^^y^'pl{a^BS+pl) 


\y±)d^i\PB, y±) 


which is the light-cone result Eq. (3.15). 

Thus, the final result for the sum of diagrams in Fig. 5 is Eq. (4.51) 


^ = -ig^P’"'' I —\{y±\^[a/3BsidiU + dlUp, 

Jrr' a 1 pi 


Pi 


, L^u^\y^r 

apBS+pl 


+ {y±\ 


uPbS 


p\{aPBS + p\) 


\yAMmB,yi.)-upy^)f] 


(4.56) 


• 2 £‘ 
= -^g f 


2 £nkl 


1 






a^BS +p]_ 


where we imposed our cutoff a > a > a'. Again, let us note that the above expression is 
valid with our accuracy in the whole range of transverse momenta. 

Similarly to Eq. (4.28) we can rewrite this formula in the form without subtractions 


T—(y±l - 4-^fc(/3s)(i +Ui) 

Ja' a Pl_ 


aPBS+p\^ 


p\{aPBS + p\) 


(4.57) 


where Fk di = diTk = —i[pi,Fk]- Indeed, in the light-cone case <C p\ one can neglect 
the operators with high collinear twist so both equations (4.56) and (4.57) reduce to the 
last terms in the r.h.s’s which are the same. Also, as we discussed above, in the shock-wave 
case {l‘]_ ~ p\) and Pb small one can replace FPPb) by Ui so the r.h.s’s of Eq. (4.56) and 
Eq. (4.57) coincide after some trivial algebra. Finally, if ~ and /3 b > ^ we have 
uPbs S> p\ so again the equations (4.56) and (4.57) reduce to the last terms in the r.h.s’s. 


4.4.5 Virtual correction for the complex conjugate amplitude 

The calculation of the virtual correction in the complex conjugate amplitude is very similar 
so we will only outline it. As in the previous Section, we start with the formula (3.16) 
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which can be rewritten as 


- / dx, = 


(4.58) 


27 r°° ^ 

^ Jxt 

1*00 ^OO . 

/ dz* / d4 -^.T(3;*,a;_L)([x*,2;*]a:A2(2:*,x_L)[2;*,4]x^2(4,a^±)K,oo]a;)”'”l 
Jx„ Jz^ ^ 


4 coo 


Using Eq. (11.28) we get 

^y dx, e-'^«"*(iC(x*,xx)[x*,(X)]r)^‘®' ^ = (4.59) 

■ 2 r°° tfry r~^* r°° _ . 

= - JLpki — / d2;*(x_L|Ue“*(^®+^)**{/3B^i,(-cx),oo;p_L) 

do ^ j — OC J CT:^ 

2 

+ -pi[^,,(-oo,cx);j?x) + Q..(-oo,oo;px)]}e*^^*|xx)^' 


+ g-ifex* / 


r°°da 

lo a 


- /3B([oo,x*]^j,(x*,z*;p_L)[^*,oo]) 


kl 


2 i 


- y{5i([oo,x*]0,,(x*,z*;4i_L)[5;*,oo]) - 5j([oo, x*])0,,(x*, z*;p_l)[ 2 ;*, oo] 

2 

- [oo,x*]0„(x*,z*;p_L)9i([2:*,oo])}^^]e“*^(*“^)*|x_L) 

f dz'^ (x_L|([oo,x*]F,i(x*)[x*,z*]0,,(2;*,2;(;p_L)K,oo])^^e"*^(^"''d.|a;_L)| 


+ 


Similarly to Eq. (4.50) it is possible to demonstrate that the last three lines in the r.h.s. 
of this equation exceed our accuracy, and moreover, one can neglect factors Using 

formulas (11.29) for and (11.10) for Q^^i, we obtain the virtual correction in the complex 
conjugate amplitude in the form 


(:Fri/3B,xx)r = ^y^x, e-*^«"*(F.7(x„xx)[x„oo]r)^ 


• Z £ 

= -W f 


2 £nkl 


"da 


I (x I 


U 


/ a I af3BS+pj_ 


{ - iafdBsdiW +Pid'j_W}^\x±) 


kl 


Pi 


+ {MldB,x±)-idiUMf'ixx\ 


a^BS 


|xx)} (4.60) 


' plia^BS + p'i) 
where we have imposed our cutoffs in a and used the formula 

/ OO 

dz* [-oo,2;*]2 

-CXO 

X (-—D^F,j{z^,z^)[z^,oo\^ + ^ / d< F,j{z^,z^)[z^,z%F,\z[,ZA_)[z[,oo\^ 

Similarly to Eq. (4.51) this expression is also valid in the light-cone case <C p\ where it 
coincides with Eq. (3.18). 
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Alternatively, one can use the expression without subtractions (cf. Eq. (4.57)) 
(-^”(/3 b,x±))" = -ig^r^^l —(xx|[/—(4.61) 


j> a apBS + 


X [25^6] - gijg^^){idk - Uk)WB)^ + Hh) 


M 


p> 


a^BS 


PI p\{aPBS + p\) 


k±) 


kl 


5 Evolution equation for gluon TMD 


Now we are in a position to assemble all leading-order contributions to the rapidity evolution 
of gluon TMD. Adding the production part (3.1) with Lipatov vertices (4.28) and (4.33) 
and the virtual parts from previous Section (4.57) and (4.61) we obtain 

\ In cr 


{J^nPB,x±)T;^iPB,y±)) 


(5.1) 


/ —Tr{/,(/3B,x^)(y^| - ^Fu{PB){idi +Ui){25t6\ - g,^g^^)U , 

J(j/ OL ^ PX. OCjjQS “r P I 


[/t 




uPbs 


p\^{aPBS+p\) 


\y±) 


+ {x±\u 


X)' ‘ 

U\26^5l - gimg^%dk - Uk)WB)^ 


aPBS + Pj_ 


Pl 


+ M^By 


a^BS 


■^\x±)Tj{PB,y±)'^ + O(a^) 


'pl(aPBS + pl}' 

where Tr is a trace in the adjoint representation. In the explicit form the evolution equation 
reads 

-AXyB,x±):Fy(PB,y±) ( 5 . 2 ) 


dlnu 

= - OsTr 

- ‘^k-tgikU 


I j d^k_i{x±\^ 


U- 


1 


aPBS+p]_ 


{ii^ku+Pktjy 


(xhsg^i - 2kj[ki 


a^BS + k?\ 


aPBS + p]_ 


— “^gfikU 


Pi 


2ky 


aPBS+pj_ 


2 + ,^2 Sik 


k] 




}X^{PB+'^)\k^) 

) as 


;/ k"^, . ( <xPbs 5^ — 2kkkj i 

X {ki_\F\PB + ^){ - / , ihU + Upi)—^ - 


crs"^l aPBS + k\ 
1 


aPBS+p]_ 




aPBS + p\ 

p 

PI 


aPBS+p]_ 
k £l ^ ^kl\ 


1 


+ 2FiPB,x^){yX - FpWB){idi +Ui){25y,5]-gjmgnU- . , 2 

a^BS+py 

uPbs 


Dt 


+d^j{^B) 


p\{aPBS+p\) 

TV 

p 


\y±) 


+ 2(xx|[7——- g^mg^%^k - iJk)WBY-2 
apBS+py Pi 


+ fpPb) 


ctPbs 


pHa^BS+pi) 


\x±)Fj{PB,y±) \ + 0{ai) 
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The operators J^j{(3) and J-i{j3) are defined as usual, see Eq. (4.29) 


{ki.\m)\yi.) = \ Idy. (5.3) 


The evolution equation (5.2) can be rewritten in the form where cancellation of IR and 
UV divergencies is evident 


d 


din a 

= - OfsTr 


(Pb, x±)J^j (/3 b , y±) 

d:‘^k±{x±\( U 


(5.4) 


al3BS + p]_ 


(U^kk+PkU^) 


(x/dBsgkLi - 2kj^ki 


- 2kj;gikU—-^^U^ - 2g^kU \ , U^}t’^{Pb + '^)\k±) 
^ a/3BS+Pj_ aons + v. J as' 


a/3BS + Pj_ 


a(3BS + k‘j_ 
as' 


1 


// k‘^, . ((x^BS^y — 2ky{kj j. , 

(k^\P /3b + ^ B + Upi y . 

us t (t/3bs + aj3BS+p\ 


,2 ^2 u^}\y±)+ 2 / d^kyxyTi{/3B + ^J\ky 


-2kk^gjlU— — -i-o ,2 

al3BS+p]_ al3BS + p]_ 

x{kx\3F\l3B + —){p - ^^^2 o \ 2 

^ as'iki a^BS + ki a^BS+pi 


k] 


as' 


+ 2U 


9jl 


-U^ - 2^U- 


Pj 


a^BS + p]_ k^ a^BS + p]_ 


+ 2 d' 


U „ ^ ^ (Uyk+PkU^)^ ^^a^^‘^^2 + 2^ 


1 


gik 


a/dBS + Pl 


2U \ u^^}t\/3b + ’^)\kx)ikx\3Fy/3B + ^)\y±) 

a/3BS+p^ k^) as' 


k\ a^BS + kj^ 

U2 


aPBS+pj_ 


■W 


+ 2JXi{/3B,x±){y±\ - ^TkiMii di +Ui)(26^S^j - gjmgyu 

1 


as' 


1 


u^\y±) 


+ 2(x±\u 


aPBS+Pj_ 


P± ' ' .. ■ (jp^g^p2^ 

17t(2hf4 - gimg^iidk - Uk)Hf3B)^\x±)Tj{/3B,y±) 


pi 


- 4 


' d‘^k_i 






+ ^,xx)J-,(/3b + 


as 


as 


a^BS 
a/dBS + k‘^. 


Ti 


(/3i3,x±)JE,(/3B,y_L) I + 0{al). 


The evolution equation (5.4) is one of the main results of this paper. It describes the rapidity 
evolution of the operator (2.8) at any Bjorken xb = /3b and any transverse momenta. 

Let us discuss the gauge invariance of this equation. The l.h.s. is gauge invariant after 
taking into account gauge link at +00 as shown in Eq. (2.7). As to the right side, it was 
obtained by calculation in the background field and promoting the background fields to 
operators in a usual way. However, we performed our calculations in a specific background 
held A, (x*, x_L ) with a hnite support in x_l and we need to address the question how can 
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we restore the r.h.s. of Eq. (5.4) in a generic field A^. It is easy to see how one can 
restore the gauge-invariant form: just add gauge link at -\-oopi or —oopi appropriately. For 
example, the terms in r.h.s. of should be replaced by Uz[z± —oopi, z'j_ — 

^Pi](-^l o-^g!j-p 2 \z')Ul,. After performing these insertions we will have the result which is 
(i) gauge invariant and (ii) coincides with Eq. (5.4) for our choice of background field. At 
this step, the background fields in the r.h.s. of Eq. (5.4) can be promoted to operators. 
However, the explicit display of these gauge links at ±oo will make the evolution equation 
much less readable so we will assume they are always in place rather than written explicitly. 

When we consider the evolution of gluon TMD (1.6) given by the matrix element (2.4) of 
the operator (2.8) we need to take into account the kinematical constraint k'j_ < a(l — I3b)s 
in the production part of the amplitude. Indeed, as we discussed in Sect. 3.3, the initial 
hadron’s momentum is p ~ p 2 so the sum of the fraction I3bP2 and the fraction -^P 2 carried 
by the emitted gluon should be smaller than p 2 - We obtain (p = Imr) ^ 


^{p\:Fn/3B,XB)T-{^B,y±)\pr 


(5.5) 


= —a 
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- 4 


„ kl 
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a^BS + ki 


^(l — /3s- —)^i{l3B H— —,x±)jFj[f3B H— —,y±)e''^’^'^ 

as as as 




{PB,x±)d'j{/3B,y±) jlp)"' + 0(af) 


' Strictly speaking, we need to consider matrix element a;x).7T(/is, J/x)|p-t CP 2 ) proportional 

to (5(5), see Eq. (2.4) 
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Note that we erased tilde from Wilson lines since we have a sum over full set of states 
and gluon operators at space-like (or light-like) intervals commute with each other. ® This 
equation describes the rapidity evolution of gluon TMD (1.6) with rapidity cutoff (2.1) in 
the whole range of 13b = xb and k± (~ |x — ylj^). In the next section we will consider 
some specihc cases. 


6 BK, DGLAP, and Sudakov limits of TMD evolution equation 


6.1 Small-x case: BK evolution of the Weizsacker-Williams distribution 

First, let us consider the evolution of Weizsacker-Williams (WW) unintegrated gluon dis¬ 
tribution (1.1) which can be obtained from Eq. (5.5) by setting 13b = 0. Moreover, in 
the small-x regime it is assumed that the energy is much higher than anything else so the 
characteristic transverse momenta ~ (x — y) 'C s and in the whole range of evolution 

(1 ;§> CT ;§> —"'^6 have ^ <C 1, hence the kinematical constraint 6(1 — Pb — ^) 

in Eq. (5.5) can be omitted. Under these assumptions, all Ti^Pb + and J-i{l3B) can 

be replaced by idiUU^ and similarly for the complex conjugate amplitude. To simplify 
algebra, it is convenient to take the production part of the kernel in the form of product of 
Lipatov vertices (4.26) and (4.31) noting that the “subtraction terms” IFi and Tj vanish in 
this limit. One obtains the rapidity evolution of the WW distribution in the form 
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din a 


un^±Miy±) 


4asTr{{x^\Up,U^{U^^m - - ^JUp,U^\yx) 


(xxlU-^^W^lxx) - ^(x_L|3^|x_L)C/i(x_L) 


Pi 


pi' 


■pi' 


Uj{y±) 


Uiix±)\{y±\^U^^U^\y±) - ^{y±\^\y±)Uj(y±)\\ (6.1) 
P_L P± ^ ^ 


where we used the formula 


-(y±l^(2<9j4C/ -FpjfcdiU) 


-u^\yA 


a (3s + p2 
p\ a(3BS+p^j_ 


o-(3bs 


■P± 


( 6 . 2 ) 


In this form Eq. (6.1) agrees with the results of Ref. [17]. To see the relation to the BK 
equation it is convenient to rewrite Eq. (6.1) as follows [24] (cf. Ref. [25]): 

^Ut{zi)U‘{z2) (6,3) 

= - + T )[/- i)T|r(c.3U - i)]{i af +f/f)} 

OTT J ^13^23 

®We have left T as a reminder of different signs in the exponents of Fourier transforms in the dehnitions 
(2.2) and (2.3). 
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where r] = In cr as usual. In this equation all indices are 2-dimensional and Tr stands 
for the trace in the adjoint representation. It is easy to see that the expression in the 
square brackets is actually the BK kernel for the double-functional integral for cross sections 
[17, 26]. Hereafter, to ensure gauge invariance, Ui{z±) must be understood as 
Ui{z±) = J^i{0,z±) = ^fdz^ [oo, oo] and gauge links at oopi must be 

inserted as discussed after Eq. (5.4). 


It is worth noting that Eq. (6.3) holds true also at small /3b up to /3b 




since 


(a:-!/)] 


one can neglect <jI3bs in comparison to 


in the whole range of evolution 1 ;g> cr 3> 

in Eq. (5.5). This effectively reduces 13b to 0 so one reproduces Eq. (6.3). 


6.2 Large transverse momenta and the light-cone limit 

Now let us discuss the case when 13b = xb ~ 1 and (x — y)^ ~ s. At the start of the 
evolution (at u ~ 1) the cutoff in in the integrals of Eq. (5.4) is ~ (x — y)^- However, 
as the evolution in rapidity (~ in cr) progresses the characteristic p\ becomes smaller due to 
the kinematical constraint p^ < a{l — I3 b)s. Due to this kinematical constraint evolution in 
cr is correlated with the evolution in p^: if cr a' the corresponding transverse momenta of 
background helds are much smaller than p^ in quantum loops. This means that during 
the evolution we are always in the light-cone case considered in Sect. 3 and therefore the 
evolution equation for 13b = x^ ~ 1 and (x — y)2^ ~ s is Eq. (3.25) which reduces to the 
system of evolution equations for gluon TMDs V{I3b, |z_l|, In ct) and 'H{I3b, |. 2 ;_l|, Inu) in the 
case of unpolarized hadron. 


6.3 Sudakov logarithms 

Finally, let us consider the evolution of'D{xB,k±,r] = In ct) in the region where xs = /3s ~ 1 

and ~ (x — p) ~ few GeV^. In this case the integrals over p‘j_ in the production part 

of the kernel (5.5) are ~ (x —?/)j^ ~ so that p^ <C (t(3bs for the whole range of evolution 
^2 „2 
1 > cr > For the same reason, the kinematical constraint 6{1 — (3b — ■^) in the last 

line of Eq. (5.5) can be omitted and we get 

J-{p\^ni3B,xx)J^-i(3B,y±)\pY^^^ (6.4) 

alncr 

= 4a,iV, [^e^^P’^-y^Hp\:F-{(3B + + ^,y±)\p) 
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As to the virtual part 
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P_L 


a/3BS + pi 


t^^|y±)|p), 


- 38 














the two last lines can be omitted. Indeed, as we saw in the end of Sect. 4.4.4, these terms 
are non-vanishing only for the region of large p\ ~ cr/Sss. In this region one can expand the 
operator O = di +Ui){2S^6^j-gjmg^^)U as 0{z±) = 0{y±) + {y-z)idiO{y±)+... 

and get 


{y±\^o 


pj_ a/3BS + pj_ 


\y±) 


Oy{y±\ 


Pm 

p\{al3BS + p\) 


\y±) + 


idm^y 

Ana^BS 


The first term in the r.h.s of this equation is obviously zero while the second is O () 
in comparison to the leading first term in the r.h.s. of Eq. (6.5) (the transverse momenta 
inside the hadron target are ~ niN ~ IGeV). 

Thus, we obtain the following rapidity evolution equation in the Sudakov region: 
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= AasNc 


{p\^i ih , ) J'" {13b , ) \p) 

d'^p± 


( 6 . 6 ) 
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ei(P’-y)x^p\:F-{i3B + + ^,y±)\p) 
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aj3BS+p]_ 


(p\^i {f^B , x±)F'^{13B,y±)\p) 


As we mentioned above, the integrals over in the production part of the kernel (6.6) 
are k‘j_ whereas in the virtual part the logarithmic integrals over are restricted from 

-| (^ •y 4 ^ 

above by an extra ^2 ^ leading to the double-log region where 1 ;:?> cr and 

ajdBS p\ {x — p)j^. In that region only the first term in the r.h.s. of Eq. (6.6) 
survives so the evolution equation reduces to 


^{p\F-{^B,xi.)F-{^B,yx)W=^'^ 

dlncT ■' 


(6.7) 


TT J pj_ ^ 


- d(p>-y)x] (pIf^Pb, x^)F^{/3b, y±)\pr 


which can be rewritten for the TMD (1.6) as 


—^T>(xs,z_L,lnc7) = - ^^^T>(a:B,z_L,ln(T) (6.8) 

a in cj TT^ JPj_ 

We see that the IR divergence at p^ —)• 0 cancels while the UV divergence in the virtual 
correction should be cut from above by the condition p^ < as following from Eq. (6.6). 
With the double-log accuracy one obtains 

—- 'D{xb, z±,lna) = - - — -V{xb, zx,l'^(y)\nasz] + ... (6-9) 

alncT TT 

where dots stand for the non-logarithmic contributions. This equation leads to the usual 
Sudakov double-log result 

V{xB,kx,\na) ~ exp { - In^ ^}T>(xb, In —) (6.10) 

^ Ztt k^ s 
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It is worth noting that the coefficient in front of In^ ^ is determined by the cusp anomalous 
dimension of two light-like Wilson lines going from point y to oopi and oop 2 directions (with 
our cutoff a < a). Indeed, if one calculates the contribution of the diagram in Fig. 6 for 
Wilson lines in the adjoint representation, one gets 



Figure 6. Cusp anomalous dimension in the leading order. 


([oopi, 0][0, 00 P 2 ]) = — ig'^NcJ dad[3d'^p± 


6{a > lal > a') 


2.. rda 
= -g^Nc — 
.L' a 


pI 


Stt^ 


Ncln^ 

a 


a{l3 — ie){af3s — p\+ ie) 
‘' dp\ 

p\ 


( 6 . 11 ) 


which coincides with the coefficient in Eq. (6.8), cf Ref. [27]. 


7 Rapidity evolution of unintegrated gluon distribution in linear approx¬ 
imation 


It is instructive to present the evolution kernel (5.5) in the linear (two-gluon) approximation. 
Since in the r.h.s. of Eq. (5.5) we already have J-k and Ti (and each of them has at least 
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one gluon) all factors U and U in the r.h.s. of Eq. (5.5) can be omitted and we get 


diner 


= — a< 


(7.1) 


iV. - fe - li) [( (7- + 

./ 1 as^ l\a/3BS + pj_ a/3BS + kj_ a0BS + p, ) 


a0BS + Pj_ 


+ ‘^dik 
+ “^gij 


kj a0BS + 2kj_ {p' + k)i 


. a0Bsd^-2k>^kj {p' + k)i _ ^k^^gp + P'j 

V a0BS + k\ a0BS+p'\ 

2gji ‘^P'jki 


PI a0BS + p 


/2 


r + 


k^ a^BS + k^ a0BS+p']_ a0BS + p'\ kj_{a0BS + p']_) 


{p + k)k kia0BS + 2k‘] 


+ 


2gik 


2pikk 


a0BS+ p‘]_k\ a0BS + k\ a^BS + p^ kj_{a0BS + p\) 2- 
2 r {2k'‘p'- - kjp'^)5^ {2pik^ - kip^)5^, i 


X {p\^^{I3b + —,P± - kx)F\0B + —,P± - k±)\p) 
as as 


+ kf (p -'+kii + 

- ^(pI ^(1 - Pb - ^)^i(PB + §7.P± - ^B)BiPB + “ Pp 


as ' 


as 
a0BS 


as 


-F 

-3 


\0B,Pi)kF'-{0B,Pi_) Ip)} 


a0BS + kj_ 

where we performed Fourier transformation to the momentum space. Also, the forward 
matrix element {p\Fi{0B,P±)J^ji0B,p'±)\p) is proportional to 6^'^\p± — p'j_)- Eliminating 
this factor and rewriting in terms of Fij (see Eq. (2.9)) we obtain {p = Infi) 


^TZij{0B,P±;g) 


= - UsNc / d' 




{2p - k)k (y0Bsgkii - 2{p - k)j^{p - k)i 


l\a0BS+p]_ 


cr0BS + (p - ky 


- 2 


(7.2) 

{p- k)j[gik+Pigij.k \ 
a0BS+p]_ 2 


a0BsS^ - 2{p - k)y{p - k)j {2p -k)i _^{P- kf^gp + (ifp. 


+ 2gik 
+ 2gij 


a0BS + (p - k)^ 
(P - k)j{2p - k)i - 2pj{p - k)i 
{p-k)\{a0BS+p\) 

{p - k)i{2p - k)k - 2pi{p - k)k 


a0BS+py 


a0BS+p]_ 


+ 


{P - k)j{2p - k)i 


+ 


2gp 


{a0BS + {p-k)\){a0BS+p\) a0BS + p\ 


+ 


{p - k)i{2p - k)k 


+ 


2gik 


X qi - & - t J 

as as 


+ 2 


- 4 


{p-k)\{a0BS+p\) {a0BS + {p-k)\){a0BS + p\) a0BS + p\ 

X 0(] 

,5f {kjp^ - 2k^pj ) + d^kip^ -2pik^) 

—7725- T^kl[PB,P±',V) 

^i±J 

a0BS 


k\[a0BS + (p - k)\] 


^(l — 0B — 


{p-k)\ 


Ipiij {0B + ,k^-,g) - 


-TZij{0B,P±]v) } 


{p-k)\ as ^ ’ k\{a0BS + k\)‘ 

Let us demonstrate that Eq. (7.2) reduces to BFKL equation in the low-x limit. Indeed, 
in this limit Fij is proportional to the WW distribution (1.1): 
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TZij{Q,k±) ~ (p|tr{i/j(x)C/j(0)}|p). In the leading-order BFKL approximation 

(cf. Ref. [14]) 


{p\iT{Ui{x)Uj{y)]\p) 

_ f d‘^q± 

47r2 J qj_ 


(7.3) 


-qiqje 




'd^q'x 

«/2 


^T{q') f 

J a 


a-^oo 27rAgg' 


Gi^iq,q') 


Here <hr(g') is the target impact factor and Gu}{q,q') is the partial wave of the forward 
reggeized glnon scattering amplitude satisfying the equation 


‘^Guj{q,q') = S^‘^\q-q')+ / d‘^pKBFKL{q,p)Guj{p, q) 


(7.4) 


with the forward BFKL kernel 

Q-sNc 


7Lbfkl(q',p) = 


TT^ 


^ -ld{q±-p±) [dp^ 


dl 


(q-p) 


Thus, in the BFKL approximation 
7Zij(0,q±;lna) = qiqjR(q±; In a) = 

and the equation for R 
d 


asqjqj f d^ 

2 vr2g2y g/2 


<hr(gO 


p'i(g-p')iJ 

duj / as \‘^ 
-ioc 27rzVgg' 


R‘^(g±;lncr) = /(i^p± ^iLBFKL(g,p)-R'"(p±;lnfT) 

diner J qj_ 


Guj{p,q') 

(7.5) 

(7.6) 


is obtained by differentiation of Eq. (7.5) with respect to Incr using Eq. (7.4). 

Now it is easy to see that our Eq. (7.2) reduces to Eq. (7.6) in the BFKL limit. As we 

k‘2 

discussed above, in this limit one may set /3 b = 0 and neglect in the argument of TZij. 
Substituting 77.jj(0, A:_l) = fei/i;ji?(A:_L) into Eq. (7.2) one obtains after some algebra 

2A:^ 


d 

diner 


i?(p_L;lncr) = 2asNc d^k± 


-?]_{?-k)l 


R{k±;lna) — 


pi 


kiip-k)^ 


•R{prMcr) 


which coincides with Eq. (7.6). We have also checked that Eq. (7.1) at p± ^ p'^ reduces 
to the non-forward BFKL equation in the low-x limit. 

Let us check now that the evolution of 


P(/3B,lner) = - ^ f d^p±TZi'if3B,P±-M(T) 


(7.7) 


reduces to DGLAP equation. As we discussed above, in the light-cone limit one can neglect 
k± in comparison to p±. Indeed, the integral over p± converges at pi ~ af^BS- On the 
other hand, extra kikj in the integral over k±_ leads to the operators of higher collinear 
twist, for example 


Jd^k^ kikj R:^{(3B,kr:lna) ~ {p\dkK{h,0±)djT'^^iPB,0±)\p)^=^^^ 

~ m^c/ij(p|J^“(/3B,0_L)J^“"(/3B,0_L)|p)'“'^ ~ m^P(^B,lncr) (7.8) 
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1^2 p ,2 2 

(where m is the mass of the target) so ^ ~ a/3^s ~ os ^ 

Neglecting k± in comparison to p_L and integrating over angles one obtains 


d 


f 


d^pi_ 7^/(/3B,P±; Incr) 


diner 


X 


J 



k‘j_{a/3BS + k‘j_) 


which coincides with DGLAP equation (3.33). 

It would be interesting to compare Eq. (7.2) to CCFM equation [28] which also ad¬ 
dresses the question of interplay of BFKL and DGLAP logarithms. 


8 Rapidity evolution of fragmentation functions 

In this section we will construct the evolution equation for fragmentation function (1.7). We 
start from Eq. (5.2) which enables us to analytically continue to negative /3s = —fip- la the 
operator form, the equation (5.2) has imaginary parts at negative /3s = —/3s corresponding 
to poles of propagators {(T/dps — ± ie)~^ but we will demonstrate now that for the 

evolution of a “fragmentation matrix element” (2.10)® 


(^“(-/3s,xs)J7(-/3s,2/±))frag = 5^(0|/“(-/?F>a;s)|p + W)(p + A|J7(-/3s,y±)|0) (8.1) 


we have the kinematical restriction cr(/3s — 1)5 > p\ in all the integrals in the production 
part of the kernel (5.2). As to virtual part of the kernel, we will see that the imaginary 
parts there assemble to yield the principle-value prescription for integrals over p^. The 


® Again, strictly speaking we should consider ^^{0\Ti{—l3F,x±)\p + X){{p + ^p 2 ) +X\Tj{—PF, 3/±)|0), 
see Eq. (2.10). 
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Tragmentation matrix element” (2.10) of Eq. (5.4) has the form 
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d\n.a 


{^i {-■, X±)F'- iV^))ba.g 
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( 8 . 2 ) 

(r^Fsg^i + 2kj^ki 
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as 


h2 
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as 
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a^FS - pi 
kl rr Pi 


t/t 


+ 2{xf\( -U 
+ 2U 
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_ 2U + 2—[/■ 

a^FS-pl kl a Pfs-pI 

.2 


f'')l!/±) 




9ik 
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as'' ' ' ' " ' as 
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■jn '' ^ 

+ 2Fii-PF,x±){y±\^M-PF)iidi +Ui){25i6\ - gj^g’^^p— ^ . 


^C^^l2/±) 


+ 2(xx|t/ 


-1 


/y-\' ' 

aU\25l5l - gimg'^%dk - Uk)Fii-M^\x±)Fj{-PF,y±) 


ajdFS — pi — ie 


Pi 


- 4 


'd^kx 
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k^, 






e{PF-l-^)M-f^F + —F±)J^i{- Pf + 


as 


as 


as 




af^FS 


+ 


af^FS 


2 ^afdFS — kl — ie a/^FS — kl + ie 


—]^i{-PF, x±)Fj{-ldF, y±) 


+ 0(ag) 


where we have restored Fie in the virtual part in accordance with Feynman rules. 

Let us prove that all non-linear terms in Eq. (8.2) can be neglected with our accuracy. 
(Naively, they were important at small (Hb but small (5f are not allowed due to kinematical 
restrictions). First, consider the “light-cone” case when the transverse momenta of fast 
fields ll are smaller than the characteristic transverse momenta in the gluon loop of slow 
fields pi ~ fci. As we discussed above, in this case with the leading-twist accuracy we can 
commute all t/’s with p_i_ operators until they form [/[/I = 1 and disappear. In this limit 
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the (8.2) turns to 
d 


2 

(-F“(-/3F,a:±)-^“(-/3F,?/±))frag = - 4asTr|(x_L|6»(/3 f “ 1 “ ^) 

g^lkp. 5'^p\p^ + 2pip'^p'^ 1 rp^(5j 6''^pj + S'jpj^ - g^jp^ 

~t~ O. //-. O.\o 0~l” 


diner 


(8.3) 


+ 2 


Pi 


(a;±l 


afdps — p\ — ie {ajdps — p\— ie)^ p\ (^Idps — p\+ ie 
9 mP\p^ + ‘^pjlPjP^ 


a/dps 


p]_{a/3ps - p]_ - ie) 


{a(dps -p\ + ie)2. 

\x±) + (y±|- 


\y±){^k(^ - Pp,X±)Ti[^ - (dp,y±)){riig 


as 


afdps 


-\y±) 


'p]_{aPps - p\ +ie)' 

X {Ti{-(dp,x±)Tj{-ldp,y±)){j.a.g^ 

Now consider the shock-wave case when ~ There are two “subcases”: when 

/dpa^: > 1 and when /dpa^ <C 1 (where u* ~ ^). In the former case we have > 1 so 

afdps ^ p\ and only two last lines in Eq. (8.2) survive. Moreover, in this case Eq. (8.3) 

also reduces to the last two lines so Eq. (8.2) is equivalent to Eq. (8.3) in this case as well. 

2 

If jdpa^f <C 1, as we discussed above, one can replace iFj{—[dp) (and J^j( — [dp + ^)) 
by Uj. We will prove now that after such replacement the r.h.s. of Eq. (8.2) vanishes, and 
so does the r.h.s of Eq. (8.3), and therefore Eqs. (8.2) and (8.3) are equivalent in the case 
of large (dpa^ also. 

Let us now prove that if we replace all Fj{—fdp) and [Fj(^ — Pp + by Uj the r.h.s. 
of Eq. (8.2) vanishes. Indeed, a typical term in Feynman part of the amplitude vanishes: 

{p + X\Ujiz)U,^ul„\0) = 0 (8.4) 

To prove this, let us consider the shift of U operator on |a*pi. Since the shift in the pi 
direction does not change the infinitely long U operator, we get 


{p + X\Uj{z)UpUl„\0) = (p + X|e*i^-“*C/j(z)17yC/]„e-*i^-“*|0) 

= + X\Ujiz)U,>ul„\0) 

which can be true only if Eq. (8.4) vanishes. It is clear that for the same reason all terms in 
the r.h.s of Eq. (8.2) (and r.h.s. of Eq. (8.3) as well) vanish. Summarizing, in all regimes 
the Eq. (8.2) can be reduced to the light-cone version (8.3). 

One can rewrite Eq. (8.3) in the form: 

'' :W(-/3F.i'±).F“(-/3F,!/±))?,il"" (8,5) 

o2 ^2 ^2 
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where we used the formula 


+ 


af^ps 


'd‘^p± 1 r afips 

2Vaj5ps — p\ +ie ' af3ps — p^j_ — ie- 


^^^9{a(3ps - p\) (8.6) 

PI 


The Eq. (8.5) is our final evolution equation for fragmentation functions valid for all (x—y)^ 
(and all f3p). 

If polarizations of fragmentation hadron are not registered we can use the parametriza- 
tion (2.10) 


zp)T-i-pp, 

= ‘^F‘^5{Cjl3pg^ - gijT^ii^F,z±,p) - A^{2ziZj + gijz\)n'l{l3p,z±,r]) (8.7) 


where PLiiPp, z±,ri) = Jd‘^k± A;_l, y) and 

'H'^{/3f,z±,p) = (Ja)" 'T-Lf^fSp, z±,p), ci. Eq. (3.26). After integration over angles similar 
to Eq. (3.27) one obtains 
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gijas'D{{l3p,z±,p) 3- + 9 ijZ±)as'H{ (Pf, z±,p) 
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/33 
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UfWf-P) /3f(^f-/3)2 15p{Pp-(3Y\ 


gijasVfiPF, z±, p) + ^{2ziZj + gijz\)asJi”{/dp, z±, p) 


where /3 = — (and a = as usual). 
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This evolution equation can be rewritten as a system (cf. Eq. (3.29) for DIS) 


A 

dr] 


OisT^ii 


1 

ZF 


z±,v) 


(8.9) 



-2 + z'{l- /) asVf^ — ,z±,r]) 
J Zp 



where z' = 1 — = 1 — Idzp- Here we introduced the standard notation zp ^ for 

the fraction of the “initial gluon momentum” carried by the hadron. By construction, this 
equation describes the evolution of fragmentation TMD at any zp and any k± ~ 

Let us demonstrate that Eq. (8.9) agrees with the DGLAP equation for fragmentation 
functions in the light-cone limit x± y±. In this limit 


diner 


asP*( —,0_L,lncrs) = —NJ I 


d 1 „ 1_X “s 


Zp 


TT 


f f dz r 1 .. ' 


( 8 . 10 ) 


z> \ -I 

X asV^ — ,0_L,lncrs) - asP^( — ,0_L,lnfJs) / -- 

Zp zp Jq 1- z') 


As explained in Eq. (3.30), with leading-log accuracy we can trade the cutoff in a for cutoff 
in . In terms of the standard definition of fragmentation functions [3] 


d^g{zF,\n p) 


( 8 . 11 ) 


we have in the leading log approximation 


df(zir,ln^2) = ^P^( —,z_L = 0,lncrs) -h 0{as) 

2i\c Zp 


so we can rewrite Eq. (8.10) in the form 
d 


as{n)dl{zF,In ip^) = 
Osifi) dz' 


dln/r^ 9 


( 8 . 12 ) 


(8.13) 


TT 


^ ([f^]+ + 7 + ^'(l-^')-2)a,(/r)d^(^,ln^2) 


easily recognizable as the DGLAP equation for fragmentation functions [22]. (Here again 
the term proportional to /3-function is absent since F9F°‘'^ is defined with an extra cts.) 
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Finally, let us describe what happens if Zi? <C 1 and we evolve from cr ~ 1 to cr 


Zp 


With double-log accuracy we have an equation 

dr] ^ Zp vr^ ^zp JPp 

(cf. Eq. (6.8)) with the solution of the Sudakov type 

, zj_, In cr) ~ exp { — In^—, 2 ;j_, In ^^) (8.15) 

Zp ' ^ 2 p Zp ^zp zj_s 

The evolution with the single-log accuracy should be determined from the full system (8.9). 


9 Conclusions 

We have described the rapidity evolution of gluon TMD (1.6) with Wilson lines going 
to -|-oo in the whole range of Bjorken xp and the whole range of transverse momentum 
k±. It should be emphasized that with our definition of rapidity cutoff (2.1) the leading- 
order matrix elements of TMD operators are UV-finite so the rapidity evolution is the 
only evolution and it describes all the dynamics of gluon TMDs (1.6) in the leading-log 
approximation. 

The evolution equation for the gluon TMD (1.6) with rapidity cutoff (2.1) is given by 

(5.5) and, in general, is non-linear. Nevertheless, for some specific cases the equation (5.5) 
linearizes. For example, let us consider the case when xp ~ 1. If in addition ~ s, the 

771^ 

non-linearity can be neglected for the whole range of evolution 1 S> cr ;§> and we get 
the DGLAP-type system of equations (3.29). If kp is small (~ few GeV) the evolution is 
linear and leads to usual Sudakov factors (6.10). If we consider now the intermediate case 
xs ~ 1 and s ^ k\^ the evolution at 1 S> cr ^ ^ will be Sudakov-type (see Eq. 

(6.6) ) but the evolution at — S> cr S> will be described by the full master equation 
(5.5). 

For low-x region k± ~ few GeV and ~ ^ we get the non-linear evolution described 
by the BK-type equation (6.3). If we now keep ~ few GeV and take the intermediate 
1 xp = Pp ^ — we get a mixture of linear and non-linear evolutions. If one evolves a 

k2 

(■f-7> rapidity) from 1 to ^ first there will be Sudakov-type double-log evolution (6.8) from 
^2 ^2 

cr = 1 to cr = then the transitional region at cr ~ and after that the non-linear 
^2 ^2 

evolution (6.3) at ^ a ^ The transition between the linear evolution (6.8) and 

the non-linear one (6.3) should be described by the full equation (5.5). 

2 

Another interesting case is xs ~ and s 3> In this case, if we evolve a 

from 1 to first we have the BK evolution (6.3) up to cr ~ ^ and then for the evolution 

m? 

between <7 ~ ^ and cr ~ we need the Eq. (5.5) in full. 

In conclusion, let us again emphasize that the evolution of the fragmentation TMDs 
(2.10) is always linear and the corresponding equation (8.8) describes both the DGLAP 
region ~ s and Sudakov region ~ few GeV^. 
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As an outlook, it would be very interesting to obtain the NLO correction to the evolu¬ 
tion equation (5.5). The NLO corrections to the BFKL [29] and BK [14, 30, 31] equation are 
available but they suffer from the well-known problem that they lead to negative cross sec¬ 
tions. This difficulty can be overcome by the “collinear resummation” of double-logarithmic 
contributions for the BFKL [32] and BK [33] equations and we hope that our Eq. (5.5) and 
especially its future NLO version will help to solve the problem of negative cross sections 
of NLO amplitudes at high energies. 

The authors are grateful to G.A. Chirilli, J.C. Collins, Yu. Kovchegov, A. Prokudin, 
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ported by contract DE-AC05-06OR23177 under which the Jefferson Science Associates, 
LLC operate the Thomas Jefferson National Accelerator Facility, and by the grant DE- 
FG02-97ER41028. 


10 Appendix A: light-cone expansion of propagators 

In this section we consider the case when the transverse momenta of background fast fields 
^_L are much smaller than the characteristic transverse momenta p± of “quantum” slow 
gluons. As we discussed in Sect. 2, in this case fast fields do not necessarily shrink to 
a shock wave and one should use the light-cone expansion of propagators instead. The 
parameter of expansion is the twist of the operator and we will expand up to operators 
of leading collinear twist two. Such operators are built of two gluon operators ~ F,iF,j 
or quark ones tp jjjitp and gauge links. To get coefficients in front of these operators it 
is sufficient to consider the external gluon field of the type A,{z^,z±) and quark fields 
x_l) with all other components being zero. 

For simplicity, let us again start with the expansion of a scalar propagator. 


10.1 Scalar propagator 

10.1.1 Feynman propagator for a scalar particle in the background gluon field 

For simplicity we will first perform the calculation for “scalar propagator” {x\pP^\y). As 
we mentioned above, we assume that the only nonzero component of the external field is 
A, and it does not depend on z, so the operator a = i-^ commutes with all background 
fields. The propagator in the external field A,{z^, z_l) has the form 


I 1 I ^ \ da 


( 10 . 1 ) 


X e 


—ia.{x—y) 


•(x_L|Pexp{ - if dz^ [— - —A,{z^)]}\y±) 
/ n I G~S S 

d 1 /* 


^°The z, dependence of the external fields can be omitted since due to the rapidity ordering a’s of the 
fast fields are much less than a’s of the slow ones. 
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The Pexp in the r.h.s. of Eq. (10.1) can be transformed to 


(x_L|Pexp{ - i / dz^[— -—A,{z^)]}\y±) 
Jy OtS S 


= (xx|e-*-(^*-*'*)Pexp{^ r*dz, (10.2) 


Since the longitudinal distances inside the shock wave are small we can expand 

F.i} - 4, {p\ DjF.,}} + ... 

as Za^s^ 

= A. - - iD^F.i) (pV - ip>D^)DjF.i + ... (10.3) 

as 

This is effectively expansion around the light ray y± + ^y*pi with the parameter of the 
expansion ~ |4| 'Cl. As we mentioned, we will expand up to the operator (s) with twist 
two. 

We obtain 


0{x^,y^;p±) = Pexpj^/ dz^ e'''^*^A,{z^)e 


(10.4) 


= 1 + — 
s 


y* 

2*5 Cj.. (2-2/)*/n ir- -T^iT- \ 


^ rA [.4. 

^ J y* 


y 

V r j.. r J../ r 4 2(z - y)* . . p . 2{z' - y). 


as 


{2p^F,i - iD^F,i) - 2- 


a^s^ 


(pV - ip>D^)DjF,, 


- ^ dz^ dz'^ [A, - 


'y* Jy* 


as 


h^F,,] [A. - 


as 


yT.,1 + ... 


It is clear that the terms ~ A, will combine to form gauge links so the r.h.s. of the above 
equation will turn to 

C)(x*,y*;p_L) = {z - y)^[2pP[x^,z^]F,j{z^) - i[x^,z^]D^ F,j{z^) 


as 
+ 2 


'y* 
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2 /•X* 


as 

rz* 


{p>p [x^,z^] - ip [x^:,z^]D^)DkF,j ) [z*, y*] 


as-" 


'y* ''y* 


H-dz* / dz^ (z - y)* z[x*,z*]F,j(z*)[z*,zjF,-^(zJ (10.5) 


- 2p>p^ ^^ ^ [x^, z^\F,j{z^)[z^, z'^]F,k{z'S)[z'^,y*] + 


as 


where dots stand for the higher twists. 

Thus, the final expansion of the propagator (10.1) near the light cone y_L + |y*pi takes 
the form 


. I 1 I , 1 •/)/ N f°^iot /■° da 

(a^ P 2 , ■ \y) = -i 0 {x*-y*) - 7 ^+i 0 {y*-xA — 

P^ + ie 1 do 2a J_^2a\ 


( 10 . 6 ) 


X e 


—ia.{x—y) 


’ix±\e 




0{x^,y*-,p±)\y±) 


Note that the transverse arguments of all background helds in Eq. (10.6) are effectively y±. 
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10.1.2 Scalar propagator for the complex conjugate amplitude 

For calculations of the complex conjugate amplitude we need also the propagator 
, , 1 


ip2_7gl2^) = 


, /■“ da 
- -y*) I 


0 


2a 


2a - 


(10.7) 


X e 


—ia{x— 


^)*(x_L|Pexp{ -i f dz^[— - —A,{z^:)]}\y±) 
Jy as s 


For the calculation of the square of Lipatov vertex we need to consider point x inside the 
shock wave and point y outside. In this case one should rewrite Eq. (10.2) as follows 

(x_L|Pexp{ - i / dz^[— -—A,{z^)]}\y±) 

Jy* ® 

= (xx|Pexp{^ rdz, (10.8) 


The light-cone expansion around x± + is given by Eq. (10.3) with y* —)■ x* 




E±I 


(10.9) 


= A,- ^-A^(2F,jp> + i&F,,) - + i&DjF,ip>) + ... 


as 


(the only difference with the expansion (10.3) is that we should put the operators yp to the 
right) and therefore 

•2^ 
s „ 

( 2 ;* — x^Y 


0{x,,y,;p±) = Pexp{^ = 1 


2^5 r*, r J 
+ — / (iz=c \_A^ — 




{2F.jp> + W^F.j) - 2- 


v2c2 


-{DjF,ipy + iDWjF,y)] 


- ^ f'dz. Fdz: [i. - [i. - yyF.,p>] + ... (10.10) 


'y* -'y* 


as 


as 


which turns to 


2ig 


0{x^,y^;p±) = [x*,y*]H- ^ dz^ [x^, z^]{2F,j{zyz^,y^y + iD^F,j{zyz^,y^] 


as 

+ 2^^^—-^[DkF,j{zyz^,y^yp’" + iDWkF,j{zyz^,yy^)}{z - x) 
8^2 

+ 

n/ , 

'y* Jy* 

(z' - x) 


as-^ 
- 2 - 


as 


dz* / dz'^ {z - x)*[x*,z*\i -iF,j{z*)[z*,z'^]Fj{z'yz'^,y*\ 

Jy* ^ 

hj{z*)[z*,z'^]F,k{z'yz'^,y*]pjp^\ + ... ( 10 . 11 ) 


and we get 


1 


da 


da 


{x\ p 2 _ ,^ y) = i0iy*-x*)l —-i9{x*-y*)l — 


le 


2a 


2a - 


( 10 . 12 ) 


X (xx \d{x*,y*, px)e-'- | yx) 

Here (in Eq. (10.11)) the transverse arguments of all background fields are effectively Xx- 
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10.1.3 The emission vertex 

For the calculation of Lipatov vertex we need the propagator in mixed representation 

1 9 

(fc|in the limit —)■ 0 where k = api + + k±: 

+ (10-13) 

First, we perform the trivial integrations over x, and x_l: 




= / dx*d^x_L e 




(10.14) 

- -d‘i]e{x - 


^9x* 


as 


.fcl ... . .p1 , . 5 




= e*"^*e*^^* 


d 

dx^—9{x - y)*{k±\0{x^,y^; k)\y±) 




-0{oo,y^,y±-,k) 

where 0(oo, y*, = G'^^^’y'^^{k±\0{oo,y^;k)\y±). In the explicit form 


0{oo,y^;y±;k) 

2ig f°° 

= [oo,y*]y-^/ dz^ {z - y)^[2k^[oo,z^]yF,j{z^,yA_) - i[oo,z^]yD^ F,j{z^,y^) 


as^ 


'y* 


+ 2— —^(A:-^A:'[oo,z*]y - ik^[oo, z^jyB^')DiF,j(z^, y±)) [z^, y^jy 


+ ^ / dz^ I *dz[ {z' - y)M[oo,z^]yF,j{z^,yi_)[z^,z%F,^{z[,yA_) 

CX.S ./ 1 /. .1 ti. '• 


as 

,2 /*oo pz* 


'y* 'jy 


- 2/c^fc^ ^ [oo, z^]yF,j{z^,yj_)[z^, z%F,i{z'^,yj_)j[z'^,y*]y (10.15) 


where the transverse arguments of all fields are y± and is replaced by k^. 
Similarly, for the complex conjugate amplitude we get 

lim /c^(x|—Tw-- \k) = e“*^^0(x*,oo,xx; A) 

fc2_5.0 — It 


(10.16) 


where 0(x*, oo, x^; A) = {x±\0{x^^co,p_\_)\k±)e or, in the explicit form 

(!)(x*, oo, xx; k) 

2ig f°^ ~ . . . _ 

= [x*,oo]a;H- 2 / [x^,z^]x{2F,j{z^,xx)[z^,(X)]xk^+iD^F,j{z^,xx)[z^,oo\x 

OiS 

+ 2 ^^~''^* {biF.,{z,,xx)[z,, oo\xk^k'- + ibWiF,j{z^,xx)[z^,oo\xb)]{z - x) 


as 

Zl roo POO 


+ ^ / dz^ I dz'^ {z - x)^[x^,z^]x( - iF,j{z^,xx)[z^,z[]xFj{z[,xx)[z[,oo\3 


OS'! 


- 2— - —F,j{z^,xx)[z*, z'^]xF,i{z'^,xx)[z'^,oo\xk^kf] + .. 


as 


(10.17) 
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In the complex conjugate amplitude we expand around the light cone x± + so the 

transverse arguments of all fields in Eq. (10.17) are x_\_. Note that the second terms in the 
r.h.s. of Eqs. (10.6) and (10.7) (proportional to da) do not contribute since a > 0 for 
the emitted particle. 


10.2 Gluon propagator 


10.2.1 Gluon propagator in the background gluon field 


As we saw in the previous Section, to get the emission vertex (10.13) it is sufficient to write 
down the propagator at x* > y*. The gluon propagator in the bF gauge has the form 


i{K{x)Al{y)) = {x\ 


1 


+ 2igF + 


rJvC 


le 


(10.18) 


xt>y, 

= —I 


— ^ (x_L|Pexp{ - if dz^ [— - —A,{z^) - —F{z^)]}\y±) 

Jy as s as 


Jo 2a 


ah 

{lU 


= ^— e 

./n 2 a 


—ia(x—y)m 


X (x_L|e 


■rl /- ^ 


Pexp( 


\ig d-z, F^^^*-yA(A, + -F)iz,)e-^^^^*-yA\ \y^) 

Jy, s ^ a ' ) yu 


where powers of F are treated as usual, for example {FA,FF )= F^A^g^^F^^Fx^. The 
expansion (10.3) now looks like 




2 

{^A,g^y-\ — Fp^y)e yF = A,g^y -\— 

a ^ a 


F„ 


+ i 


.Zif y* T 2 


as 


— 9y,v 


9uu “1“ r yy 

a 

[p±,A,gf_i,y + -Fyu] - ^20,2^2* bi’ [p±,A,g^y + -Pyu]] + ••• 

-P - - ip^DkD^F.,) 


A, - — —{2]FF,j - iD^F,j) - 2 


as 


a^s 

+ -F^y + 2i^F^p)DjF^y + 2i ^^~^yf* jFp^D,DkF^y + ... 


a 


a^s 


(10.19) 


- 53 











so we get 


Gij.uix*,y*;p±) = Fexp^ig f d‘^z^e^<^sF* ?/.)2^*)| (10.20) 


+ ig Hd-z, + -F) (z*)e-* - - g'^ Hd-z, 

Jv. s ^ a s 

cz* O 


— 9^lu 


(A, + -F) Hd-zi + -F)^ (Oe"*- (^•“2'*) 

a s a " 


^ ^ r:, / 

- “1“ / I 

■5 Jj/* ^ 




M^*) - (2p^F., - iH'F.,) - 2^^-^{p>p’^D,F,k 


- ip^DkD^F.j 


+ —F^y + 2i 


as 
.z^-y^ 




a 


p _|_ 2i— _ 


a^s 


p'p’^DjDkF^y} 




'y* -Jy* 


a 


as 


X [5«^. + lFi-2/^^^-^F.fe5«] + 

a as 


+ 2 


2i 
as 

“ y)* 


f * / 2i 

= 9tMAx*,y*]+ g dz^ i - ;^{z-y)^:gf_,y{2p’[x^,z^]F,j{z^)-i[x^:,z^]D^F,j{z^) 

Jy* ^ 


as 


(p^p^[x*, z^]DjF,k - ip^[x^,z^]DkD^F,j)} 


+ -^{di,P2i' - 5lp2y){[x^,z^]F,j{z^) + — — —p^[x^, Z^]DkF,j{z^) 
0^5 G~S 


+ 2 


iz-y)lkj 


a^s^ 


P P , Z:^ ]DkDiF»j }^ [z*, y*\ 
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2 pXt yz* 


as-" 


'y* -Jy* 


H- 3 / dz^ dz'^ ( [igfj,u{z' - y)^ - P2yP2u] [x*, z^]F,j{zZ[z*, z'^]F,\z'^) 


as 


- 2^p>p^{z - y)^{z' - y)*[x^, z^]F,j{z^)[z^, z'^]F,k{z'^)'^ [z'^,y*\ 


Note that Ff^^F^"^F^y and higher terms of the expansion in powers of Ff^y vanish since the 
only non-vanishing field strength is F,i. 

Finally, 


(x| 


1 


+ 2igF + 


le 


’ da 


rJy)yiy = -id{x*-y*)l -pp + iO{y* - x^) I — (10.21) 


2a 


da 


X e 


-^-^--y)-{xx\e-"^^^-y^*Qfy{x.^yppx)\y^ 


For the complex conjugate amplitude we obtain in a similar way 


(^i „, .„L ,j i/C = [»(!/.-^.)^ 


' P2 + 2igF - 


le 


( 10 . 22 ) 


X e 


—("-^)*(xx|g“^(x„y,;px)e-*^("-^)*|yx) 
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where 


Q^J.l'{x*,y*■,p±) = Fexp^igJ [A, + ^F){z^:)e ^*)| (10.23) 

fv* / 2i ~ . _ . _ 

= gij.u[x*,yA + 9 I dz^ [x^,z^]i — ^{z - x)^gf_iu{2F,j{zA[z*,y*]p^ + iD^F»jiz*)[z*,yA 

/ /Y> ' G~S 


+ 2 


s 
2i 
as 
{z-x) 
as 


(DjKk (zAlz* , y*\p^p’" + iDk&F,j (z*) [z* , y*]p^)} 


- -^idiP2u - dip2ii){F,j{zA[z*,y*] + DkKj{z^ ) [z^ , y^P^ 

+ 2 ^^~^f* DkDiF,j(zAlz,,y,W^ 


a^s^ 


8^2 ry, f^* { 2 ~ 

H-3 / dz^ / dz'Ax*, z^]i^[- ig^t,{z - x)^ - P2^,P2u]F,j{zA[z,:, zi]F,Azi)[zi,yA 

t/ X* J Z* OiS 


- “ x)*(z' - x)*F,j(z*)[z*,2;(].F,fc(4)K)y*]F'P^) 

10.2.2 Gluon propagator in the background quark field 

We do not impose the condition D^F,i = 0 so our external field has quark sources D^F^- = 
gil>F;fii'ip which we need to take into consideration. The corresponding contribution to gluon 
propagator comes from diagrams in Fig. 7 



Figure 7 . Gluon propagator in an external quark field. 


(.4;{x)4(s)>Fig, 7 = g^liAp^\^r{Ap^^\y)'^ (10.24) 

X [{z\'^t''ik^ iP p2 1 ■ 1 ■ 

As we mentioned above, we can consider quark fields with spin projection onto pi 
direction which corresponds to 'i/i(...)V’ operators of leading collinear twist. In this approx¬ 
imation P 2 A’ = 0 so the only non-zero propagators are (A,(x)A,(y)), {A,{x)Ai{y)) and 
{Ai{x)Aj{y)). In addition, we assume that the quark fields '4>{z) depend only on z± and 
z* (same as gluon fields) so the operator a = commutes with all background-field 
operators. We get 

{A^.ix)Al{y))F,,. 7 = 2tg^ [d^zd^z' (x|-^|z)-(z'|-^|y)''^ (10.25) 

J + le + le 

X [{z\i)f^lP, ^ . + {z\i}t^ P, ^lfll2\z)] 

F le F Fie 
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In our gluon 
rewrite as 


field yP = a + 2^P,+ j5_l so = 2q;P, — p\+ i> 2 l^F,j and one can 


P, 


2aP, — p\+ie 


.2a 2a 2aP, -p\+ie 


1 


^ , A _ 

2q/ o^, o^, d /V-.2 


^u; 2a 2aP, — p\+ie 
(the term | F»j does not contribute due to = 0). Similarly. 

1 


(10.26) 


so 


2aP, — p\+ie 
one can rewrite the propagator (10.24) 


P,= — + -^ 

2a 2aP, - p\ 


Pi 


— p\ + ie 2a 


(10.27) 


as 

1 


-D^F„ 


1 


-\yy 


ab 


(7l“(x)7l^(y))Fig. 7 - ig{x\ ^^^2 ^ r,, ^2 ^ 

'“"(pi Vi + 2ip>dji> - di^yy p2\-^ ^AA p2\iX ^\y) 


(10.28) 


^a(P^ + 7 e)^ " " “ ' i^^ + ie^ ' + le' ' ' 

+ ^ a{P^+ ApI - 

where in the first line we have rewritten '^I'lp as —g{D^F,jY^. 


Similarly, we get 


(A“(x)4^(y))Fig. 7 (10.29) 

= - - ^dAhj ^I7i^" p2 p 

(7lf(x)7l5(y))Fig. 7 = ( p2 p ^i7i^" p2 p i/ '^^ p2 + iX ^\y^ 

for the remaining propagators. 

If now the point y lies inside the shock wave we can expand the gluon and quark 
propagators around the light ray y_\_ + It is easy to see that the expansion of the 

gluon fields A, given by Eq. (10.3) exceeds our twist-two accuracy so we need only expansion 
of quark fields which is 

2 2 

^ + (10.30) 

as 

(and similarly for Y and D^F,i). 

It is convenient to parametrize quark contribution in the same way as the gluon one 

( 10 . 21 ) 


da 


da 


(7l“(x)A^(?/))Fig. 7 = -0ix^-y*)l ■^+0{y*-xYl 7^ 


2a - 


^-ia(x-y). 


2 2 
t ^-L f rr> _ n ^^/-§ h\ / \ I \ . / I ^^/~i hi / \ _ 0 i 


{x±\e ^^*Qyiix*,y*\p±)\y±) + {y±\Qy!l{x*,y*;p±)e 2^)*|x_l) ( 10 . 31 ) 


flU \ 
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In the leading order we need only the first two terms of the expansion (10.30) which gives 


a^s 


= - I dz^{[x^, z^]D^F,j{z^)[z^,y^] + ^ p^[x^, z^]DkD^F,j[z^,y^]y^ 


as 


+ 


g 


2 rxt 


a^s'^ 


pXt r _ . _ 

/ dz^ dz^ {p\^l;{z^:) + 2ip>lp Dj (z^)) fil[z^:,x^:]t°'[x^:, y*]t’’[y*, zi]'ilj{zi] 
Jy* Jy, ^ 


+ 2 


+ 2 


Zif y^f 2 


as 


p\p^ ?/) Uj (z*) [z*, X*][x*, y*][y*, 4]V’(4) 


z*-y* 2 


as 


p\pPi2{z^) ]ii[z^,x^]F[x^,y^]r[y^, z'^]DjTp{z'^) 


(10.32) 


and 


Q,,{x*,y*]Px) = 

^2 rx* r^* 


9 


a^s^. 


y* 


r^* r _ 

dz^ dz'^ '^{z^)]^i[z^,y^]t^[y^,x^]F[x^,z'^](^l^{z'^)p\-2iDj'il;{z'^)p’) 
J z* *- 

V’L>j {z^) 4]V’(4)piF' 


- 2 


- 2 - 


2 * 

2* - y* 


as 


as 


-i^{z^)A[z*,y*\t''[y*,x^]F[x^,z'^]Dj'ilj{z'^)p\jP (10.33) 


Similarly, for propagator {A,{x)Ai{y)) one gets 


Q.i {x*,y*]Px) 

9^ 


a^s‘^ 


+ 


+ 


Jy, Jyt, 

as 

2iz'-y) 


f^* r _ _ ^ 

dz^ / dzi Ip’ijiz^) - iip {z*)hj ^iJi[z*,x^]F[x^,y^]t’'[y^, z'^]i){z'^) 

J V* ^ 


p^p^tpDk {z*)jj ■^i7i[z*,x^]F[x^,y^]r[y^,z'^]ip{z'J 


as 


-p>p^il)-1j ^i7i[z*, [x*, y*]t^[y*, z'^]Du'il^{z'^) 


(10.34) 


Qt{x*,yAP±) 

g2 r^* r^* 


g2 rx, j-Xt 

dz^ dz'^ V’(^*)7j^i7i[^^*,2/*]i^[y*,a:*]t“[x*,z(](V'F'+*^V)(4) 

^ ^ Jyt Jzt, '- 


- 2 


- 2 


z^ -y^ - 


as 

- y* 

as 


il^Dk {z*hi I^ilj[z*,y*]F[y*,x^]F[x^,z'^]ip{z'^)p>p’^ 


i^iz*hi hi3 [z*, y*h[y *, x*]t“[x*, zl]Dki){z'^)pPp^ 


(10.35) 


For the propagator {Ai{x)A,{y)) the corresponding expressions Qfhx*Ty*'iP±) 

Qf^(x*, y*; y_L) are obtained from Eqs. (10.34) and (10.35) by replacements in the r.h.s.’s 
Ij hli li hlj and 7 i |ii 7 j 7 ^ hn, respectively. 
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Finally, for the propagator {Ai{x)Aj{y)) we obtain 




2 l•x^, pz, _ 

/ dz^ / dz^ 

Jv* Jv» ^ 


-IJ 


as^ 


+ —— —p’^'ipDk (^*)7i^i7i[^*)^*]^“[a;*,?/*]t^[y*,4]V’(4) 


+ 


as 

2{z' -y)^ . - 


as 


P^i^li^ili[z*,x^]t°-[x^,y^]t\y^,z[]Dk'4){z[) (10.36) 


^2 rxt: 


Qij{x^,y^-,p±) = - 2 / dz^ dz'^ ipiz^hj'^Ili[z*,y*p[y*,x^]t°-[x^,z'^]i){z'^) 

J J L 


- 2 


as 

zi-y* T 


as 


'fpiz^hj A'yi[z*,y*]r[y*,x^]t'^[x^, z'^]Dk'il^{z'Jp’" (10.37) 


For the complex conjugate amplitude we get in a similar way 

r r°° d'a 7° (tai 

(W4(!/)>Flg,T =[-»(!/. -I.)/ -^+e(x.-y,)l —^ 


/O 


2a. 


^-ia{x-y). 


(a:±|Q“t(a;*, 2 /*;P±)e |y±) + (y±|e ^^*Q““(x*,y*;p_L)|x_L) (10.38) 


7l, 


• ab 


where 


Q'^lix*,y*]P±) 

rx* 




+ 


a^s 


a^s^ 


'y* 


as 


ry* 

7^* r 

/ (i2;* 

/ 

/ X* 



-ipiz^) Mz*,x*]t‘'[x*,y*W[y*, z'J{'ip{z'^)p‘i - 2iDjil;{z[)p>) 


+ 2 


2 ;* — X* 7 


as 


i^Dj (z*) ?/*]r[y*, 2;(]V’(2;')pi7^ 


+ 2 - 


. z* - X* 7 


as 


( 2 ;*) [ 2 *, X* ] t“ [x*, y* ] t'’ [?/*, ] £»j -0 (4 )piF' (10.39) 


“ ab 


Q,,{x*,y*-,px) 


a^s^ 



- 2 


2;* — X* 


as 



as 


(p^V’(^*) + 2ip>ipDj ( 2 ;*)) ^i[ 2 ;*,y*]t^[y*,x*]t“[x*, 
___ — 

p]_p>'4>Dj ( 2 ;*) ^ 1 ( 2 ;*, x*]t“[x*, 2 :(]Vi( 2 ;() 

( 2 ;*) [ 2 ;*, y*][p*, X*]t“ [x*, 2 ;(]( 2 ;() 


(10.40) 
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and 


Q‘^tix*,y*;p±) 

-.2 


-9 




fV* fV* r ~ _ . _ . 

/ dz^ / dz[ V’(2^*)7i y*]t[y*^z[]{'^p>+iD^i)) 

J J Z:^ *- 

__ i — 

4>Dk {z*)lj :^i^i[z*, x^]t°-[x^, zl]^{zl)p’ 


+ 2 - 


' Z* 

. X* T 


+ 2 - 


as 

Z' — X* T 


as 


i’{z*)lj ^i7i[z*,x^]t‘^[x^,y^]f[y^,z'jDkip{z'^)p>p’^ (10.41) 


-g 

a^s^. 


“ ab 

Q,i {x*,y*]p±) 

pi ry* rzt 

/ dz^ / dzl 

'X* J 
2 

as 

_ o < - a:.. 
as 


- ii)V)7i 7i [2:*, y* ] i^ [y*, a:* ] t] Vi (2;() 

yVi Dk {z^)ji]ii 7j[2;*, y*][y*, X*]i“[x*, ^(2;() 

yp*'?V(^*)7iA7j[^*>y*]i^[y*>a:*]t“[x*,2;(]Zifc^(z() (10.42) 


To get Q“, one should again make the replacement pj ^ 17 * —)• 7 * ^i 7 j in Eq. (10.41) and to 

“ ab 

get Qj, the replacement Pi'^iPj —)• Pj'^iPi in Eq. (10.42). Finally, similarly to Eq. (10.36) 
one obtains 


Qijix*,y*-,p±) = I dz^ I dz'^ 'ip{z^)pi^i'yj[z^,x^]t°'[x^,y^]t^[y^,z'^]i) 

'X* J Z:^ *- 

, 2^* X* — 


2 ry* 


ry* 


-'ll v*^=t=; y=t=5/^-L/ 9 

as^ 


+ 2 - 


'-ipDk (^*)7i A7i[^*)2:*]t“[x*,y*]t^[y*,2;(]V'(4)y^ 


+ 2 - 


as 

-Z'-X. 


as 


'ip{z^)pi ^i7j[2;*, x*]t“[x*, 2;(]T)fcV'(4)y^ (10.43) 


= ab ry* 

Qii{x*,y*]P±) = - n dz^ 


dz' 


- 2 


as 

2;* — X* 


V' 7 i 7 i [2*, y* ] i^ [y*, X* ] t“ [x*, ] ^ ( 2 ;() 


as 
z'-x 


p^i) Dk {z^)pj '^iPi[z^,y*\r[y*,x^]t°‘[x^,z'^]'ip{z'^] 


- 2— - -p ^(z*)7j I^i7ilz*,y*]i ly*,x*jt''[x^, z'jDki/-’(z'J (10.44) 

as J 


10.2.3 Final form of the gluon propagator 

Assembling terms from two previous Sections we get the final result for 
background-Feynman gluon propagator in external field in the form 


(A^^(x)At(y)} 


- y(x* - y*) 


''°°da d:a 

+ y(y* - X*) i 




2 a 


2 a - 


,-ioi{x-y).^ 


e "'"{(xxle "as 


cts ' 


(10.45) 


X [Qfy{x*,y*]Px) + QfAx*,y*;pi_)]\yx) + {yx\Qfu{x*,y*]Pi_)e ^^*|x_l)} 
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for Feynman propagator and 


{Ki^)Aliy)) 


(10.46) 


- 0{y^, - X*) 


'‘°°da da 

+ 0{x^: -y*) I 


10 


2 a 


2 a - 




for the anti-Feynman propagator in the complex conjugate amplitude. 


10.3 Vertex of gluon emission 

Repeating the steps which lead us to Eq. (10.14) we obtain 

lira k‘^{A'^{k)Al{y)) = - id^yO’l,l{oo,y^,y±] k), (10.47) 

01,lioo,y^,y±-,k) = G'^lioo,y^,y±]k) + Q'ji{oo,y^,y±]k) + Qli{oo,y^,y±;k) 

where 

Gfu{°o,y*,y±]k) = d''’‘’y'>^{k±\g'^l{oo,y^-,p±)\y±), 

Q‘li{oo,y^,y±;k) = d^’^’y^^{k±\Qf^{oo,y^,p±)\y±), 

Qfuioo,y.,yr,k) = e-*("’^)^(2/^|Q“t(oo,y*;p±)|fe^) (10.48) 

The explicit expressions can be read from Eqs. (10.20) and (10.32) - (10.37) by taking the 
transverse arguments of all fields to be y± and replacing the operators p> with k^ similarly 
to Eq. (10.15). 

Similarly, for the complex conjugate amplitude the emission vertex takes the form 

lim /c^(^“(x)Jiy(A;)) = oo, x_l;/ c), (10.49) 

>0 ^ ^ 

0“^(x*,oo,x_l;/ c) = G'^l{x^,oo,x±]k) + Q1;ll{x^,oo,x±]k) + Q^^{x^,oo,x±;k) 

where 

G'^l{x*,oo,x±]k) = e~^^’‘’^'>x{x±\g'^l{x^,oo;p±)\k±), 

Q‘li{x*,oo,x±]k) = e“*(*^’"')-L(x_L|Q“t(x*,(X);p_L)|A;±), 

Q^^{x^,oo,x±;k) = d^’'’''^^{k±\Q^^{x^,oo;p±)\x±) (10.50) 

Again, the explicit expressions can be read from Eqs. (10.23) and (10.39) - (10.44) by 
taking the transverse arguments of all fields to be x_l (and replacing the operators p^ with 
y) similarly to Eq. (10.17). 


11 Appendix B: Propagators in the shock-wave background 

In this section we consider propagators of slow fields in the background of fast fields in the 
case when the characteristic transverse momenta of fast fields (A:_l) and slow fields (Z_i_) are 
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comparable. In this case the usual rescaling of Ref. [8] applies and we can again consider 
the external fields of the type with Ai = = 0. 

Actually, since the typical longitudinal size of fast fields is u* ~ ^ and the typical 

A 

distances traveled by slow gluons are ~ ^ our formulas will remain correct if /i > k]_ 
since the shock wave is even thinner in this case. As we discussed above, we assume that 
the support of the shock wave is thin but not infinitely thin. For our calculations we need 
gluon propagators with both points outside the shock wave and propagator with one point 
inside and one outside. It is convenient to start from the latter case since all the necessary 
formulas can be deduced from the light-cone expansion discussed in the previous Section. 
To illustrate this, let us again for simplicity consider scalar propagator. 


11.1 Propagators with one point in the shock wave 

11.1.1 Scalar propagator 

For simplicity we will again perform at first the calculation for “scalar propagator” 
{x\-pT^\y)- As usual, we assume that the only nonzero component of the external field is 
A, and it does not depend on z, so the operator a = commutes with all background 
fields. The propagator in the external field A,{z^,,z±) is given by Eq. (10.1) and (10.2) 
which can be rewritten as 


{x\ 


1 


P2 + ie 


7:\y) = 


X 


-ie{x^-y^)f ^+i9{y^-x^)[ ^ 
Jo 2q; ./_,^2aJ 


' —CXO 


^-ia{x-y). 


{x±\e 


Pexpj^^y dz^e^cJsi^* ^*^A,(z*)e yA 


( 11 . 1 ) 

}l2/±) 


Suppose the point y lies inside the shock wave (the point x may be inside or outside of the 
shock wave). Since the longitudinal distances 2 ;* inside the shock wave are small (~ ^) 

we can use the expansion (10.3) but the parameter of the expansion is now ~ ^ <C 1 
rather than twist of the operator. Consequently, the last term in Eq. (10.3) can be neglected 

p2 

since it has an extra factor —it* in comparison to the second term: 


e^^A*-y*)A,e-^^A*-yA 

= A.- ^AJZyi(^2p^F,i - - 2. ^^* ~ - ipPD^)D,F,i + ... 

as a^s^ 

= A.-^A^(^2p^F^._iDiF,i) + ... ( 11 . 2 ) 

as 

This is again the expansion around the light ray y± + ^y*pi but now with the parameter 
of the expansion ~ ^ 1- However, we need to keep the second term of this expansion 

since the first term forms gauge links (for example, it is absent in the A, = 0 gauge). 

Since there are no new terms in the expansion (11.2) in comparison to (10.3) we can 
look at the final result (10.5) for 0(x*,y*;p_L) and drop the terms which are small with 
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respect to our new power counting. This way the Eq. (10.5) reduces to 


0{x^,y^;p±) 


= [x*,y*] - 2 / {z-y)*{ 2 p’[x*,z^]F,j{z^)-i[x^,z^]D^F,j{z^)}[z^,y^] 

as Jy^ 

+ -^ [ d^* [ dzi {z -y)4x*,z^]F,j{z^)[z^,z^]Fj{zi)[zi,y^] +... (11.3) 

J V* j V* 


'y* -'y* 

and the propagator has the form (10.6) 


(^1 p2^ • 1^) = \-i0{x^-y*)[ ^ + id{y*-x^)[ ^ 

+ ie L Jq 2a J_^ 2a J 


(11.4) 


X e 


P 


0{x^,y^-p±)\y±) 


As we mentioned, this formula is correct for the point y inside the shock wave and the point 
x inside or outside. 

Similarly, for the complex conjugate amplitude we obtain the propagator in the form 
(10.12) with 


0 {x*,y*]p±) 


= [x*,y*] + I dz^ (z-x)^ [x^,z^]{ 2 F,j{z^)[z^,y^]p>+i&F,j{z^)[z^,y^]} 


ry* 


8i9 

as 


as 

2 r^t r^* 


r^* 

^ dz^ dz'^ {z - x)4x*,z^]F,j{z^)[z^,z'^]F,^{zl)[z'^,y^] +... (11.5) 

Jv* J Vt 


which is the expansion (10.11) but with fewer number of terms. Again, the formula (10.12) 
with 0(x*,?/*; p_l) given by the above expression is correct for the point x inside the shock 
wave and the point y inside or outside. 

The expressions for particle production are the same as (10.14) and (10.16) with 
0{oo,y^:,y±] k) and 0{x^,oo,x±]k) changed to Eqs. (11.3) and (11.5), respectively. 


11.1.2 Gluon propagator and vertex of gluon emission 

As we saw in previous Section, the gluon propagator with one point in the shock wave can 
be obtained in the same way as the propagator near the light cone, only the parameter of 
the expansion is different: rather than the twist of the operator. Careful inspection 

of the expansions (10.19) and (10.30) reveals that there is no leading or next-to-leading 
terms with twist larger than four so we can recycle the final formulas (10.45) and (10.46) 

p2 

for gluon propagators. At (^<7* <C 1 the expression (10.20) for Q^u{x*,y*]Pi_) turns to 

G^lu{x*,y*]p±) = (11.6) 

/ ‘2i 

= gfj,u[x*,y*] + g dz^ i - ^{z - y)^gy,u{ 2 p’[x*, z^]F,j{z^) - i[x^, z^]D^F,j{z^)} 

Jy* ^ 

+ ~^(dl,P 2 u - 5ip2f^)[x^, z^]F,j{z^)^[z*,y*] 

gg2 rXi, j-Zt 2 

+ dz^ dz'^ [ig^v{z'-y)^ - P2iiP2v\[x*,z^]F,j{z^)[z^,z'^]F,\z'^)\z'^,y^] 

as-^jy, Jy, as 
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Looking at quark formulas (10.32) - (10.37) we see that at ^<7* <C 1 the only surviving 

terms are the first terms in the r.h.s’s of these equations. Let us compare now the size of 

these terms to the gluon contribution (11.6). The “power counting” for external quark fields 

in comparison to gluon ones is dz^ijj ~ dz^D'^F,i{z^) ~ I'jjJ ~ l\ and each 

extra integration inside the shock wave brings extra cr*. The first lines in r.h.s.’s of Eqs. 

2 / ^ 

(10.34) and (10.35) are of order of Jdz^ippi'il:{z^) ~ so they can be neglected 

in comparison to the corresponding term ^fdz^ F,i{z^) ~ in Eq. (11.6). As to the 
terms (10.36) and (10.37), they are of the same odrer of magnitude as next-to-leading terms 
~ in Eq. (11.6) so we keep them for now. With these approximations we obtain 

Qliix*,yFP±) = -^^P2 imP2u[ dz^{[x^,z^]D^F,j{z^)[z^,y^])'^’’ (11.7) 

as Jy^ 

g2 rx, r-z, 

Jy, Jy, 


Qfuix*,yFP±) = - 


-'fiiy 


as^ 




and the gluon propagator is given by Eq. (10.45) with the above and Qyu'- 

MJWAtfa)) (11.8) 


'' d:a , , f (fai ^_ia{x-y). 


-9{x^-y^) - - h6'(y*-x*)/ — 


2a 


2a. 


e 

2 


X [Gyt{x*,y*;p±) + Qfuix*,y*;p±)]\y±) + {y±\Qyl{x^,y*-,P±)e ^^*|x_l)} 

As in the scalar case, it is easy to see that Eq. (11.8) holds true if the point y is inside the 
shock wave and the point x anywhere. 

Similarly, in the complex conjugate amplitude the gluon propagator is given by Eq. 
(10.46) with 

Qy,v{x*,y*]px) (11-9) 

ry* / 2i ~ . _ . _ 

= gyy\x*-,y*\F g \ dz^ — i^{z-x)^gy,u{2F,j{z^)[z^,y^]p>+iD^F,j{z^)[z^,y^]] 

J X, \as 

4 . . 

- 2 - dip2y)F,j{z^)[z^,y^ 


as 


2 M 

F ry 


+ ^ [ *dz^ [ *dzl[x^, z^]([- igy^{z - x)^ - —p 2 ^,P 2 Iz]F,j{z^)[z^,zl]FJ{zl)[zi,y^]] 

O^S .f 'Y' ./ y O^S / 


and 


Qfuix*^y*'^P±) = ^^P2iiP2u I dz^{[x^,z^]D^F,j{z^)[z^,y^]) 


ry* 


ab 


a^s 
2 ry. 


( 11 . 10 ) 


+ ^/ I dz'^i^{z*)-fy ^i'y^[z*,x*]F[x*,y*]t^[y*,z'^]'ip{z'^), 


ry* 


as 

z ab n2 ry* 


Qnu{x*,y*-,px) = dz* dz'^ ^i*yy[z*,y*]f[y*,x*]F[x*,z'^]'ij){z'^) 


as^ 
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The expressions (11.9) and (11.10) are valid for point x inside the shock wave (and point y 
inside or outside). 

The corresponding expressions for the Lipatov vertex of gluon production are given by 
Eqs. (10.47) -(10.50) with Q, Q, and Q changed accordingly. 


11.2 Propagators with both points outside the shock wave 

In this section we will find the propagators with both points outside the shock wave. Again, 
we assume that the characteristic shock-wave transverse momenta are of order of transverse 
momenta of “quantum” fields with a > a'. As discussed in Sect. 2, we consider the width 
of the shock-wave to be small but hnite, consequently we can not recycle formulas from 
Ref. [8] for the infinitely thin shock-wave. 


11.2.1 Scalar propagator 

As in the previous Section, for simplicity we start with the scalar propagator (10.1) 


(x| 


^ LO _ • f ^ -ia{x-y) 

+ Vo "" 


2 a 


•(x_L|Pexp{ -i f dz^[— - —A,(z*)] }|y_L) 
Jy, ^ 

( 11 . 11 ) 


The Pexp in the r.h.s. of Eq. (11.11) can be transformed to 


(x_L|e ^ o^^*Pex.p^ig J d-z^ A,{z^)e * ai |e* |y_|_) = J d'^z±d‘^z'j_ (11.12) 


d. 


X (x_L|e <^‘>^*\z±_ 


_)( 2 ;_L|Pexp|ic/ j 




d-z^ <=‘P*A,{z^:)e 


A 




A, 


'\y±) 


Next, we use the expansion (11.2) at = 0 


2 2 

= A, - —{2p^F,i-i&F,i) + ... 
as 


(11.13) 


This is an expansion around the light cone z_\_ + -z^pi with the parameter of the expansion 

p2 

~ -^cr* <C 1. Note that similarly to Eq. (11.2) we need to keep the second term of this 
expansion since the first term forms gauge links. 

From Eqs. (10.4), and (10.5) we obtain (cf. Eq. (10.6)) 


, I 1 I . 1 ././ \ f°° -nr \ V 

{x\ p 2 I . \y) = -i0{x^-y*)l tt + - x^) I — 

F^ + ie 1 Jq 2a 2a J 


^-ia{x-y). 


(11.14) 






X d z±{x±\e "<-’>''*0{x^,y^;p±)\z±){z±\F‘-sy*\yj_) 


where 


0{x^,y^]p±) 

= f dz^ z^[2p>[x^,z^]F,j{z^) - i[x^, z^]D^F,j{z^))[z^,y^] 

as Jy^ 

+ [ *dz^ [ *dz'^ z'^[x*, z^]F,j{z^)[z^, z'^]F,^{z'^)[z'^,y^] + ... (11.15) 

Jyt, J y* 
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Here the transverse arguments of all fields turn effectively to z±. Note that this expression 
is equal to Eq. (10.5) at = 0. For the complex conjugate amplitude one obtains (cf. Eq. 
( 10 . 12 )) 


(^1 p 2 ^ ■ 1 ^) = ^- i0 { x ^- y ^)[ ^ 

-le 1 Jq 2a 2a J 


^-ia(x-y). 


(11.16) 






X I d z^{xi_\e *“-"'*0(x*,y*;p_L)|2;_L)(2:_L|e*-^^*|y_L) 


where 

0{x^,y^-,p±) = - ^ dz^ z^[x^, z^]{2F,j{z^)[z^,y^]p> + iD^F,j{z^)[z^,y^]) 

Oi,S .} -T: 


8ig 


,•^2 ry* rv 


as-^ 


* fy* 

dz^ / dz'^ z^[x^, z^]F,j{z^)[z^, z'^]:,FJ{ z'^)[z'^,y^] 

J Z* 


(11.17) 

f J 

Again, this expression can be obtained from Eq. (10.11) by taking x* —>■ 0 in parentheses. 

11.2.2 The emission vertex 

Similarly to Eq. (10.14) we get 


- \y) = / dx*dV 

fc2^o + ie J 


(11.18) 


■fl 

i _i. ' 




X - -^dljeix - y)^ [d^z^ (x_L|e *°i^*|z±)(^;±|0(x*,y*;p_L)e*-t^*|y_L)e*"^* 


/ o 2 

dx^—9{x - y)*{k_i_\0{x^,y^;p±)F^y*\y±)F‘^y' 


= Jd^zs. e-*(^’^)^0((X),2/*;2x;A:)(zx|e*-^*|yx)e“^* 

where 0{oo, y*,z±] k) is obtained from Eq. (11.15) 

0{oo,y^,zr,k) = F^^’'''>^{ki_\0{oo,y^]pi_)\zi_) 


(11.19) 


2ig 

= [oo,y*]^-2/ dz^ z^[oo,z^]z{2k^F,j - F,j){z^,z±)[z^,y^], 

as Jy^ 


+ 


8ig 


as-^ 


2 roo t-Zt 

/ dz^ dz'^ z'^[oo,z^]zF,j{z^,z±)[z^,z'^]zF,^{z'^,z±)[z'^,y^], 
'y* J y* 


For the complex conjugate amplitude we get 

lim A;2 (x|^^^‘^- \k) = e"* 

fc2^>o P^ — ie 


{ 11 . 20 ) 


= J dy^-^6{y - x)^{x±\e~"^'"*d{x^,y^;p±)\k±)e~"°‘'"‘ 
= Jd^zx e*(^’")^(xx|e-*-^*|zx)0(x*,(X);zx;A;)e-*“^* 
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where 0{oo, y*,z±-, k) is obtained from Eq. (11.17) in a usual way 
oo; z_l;/ c) = e“**^^’^^^( 2 ;_L| 0 (x*,oo;p_L)|/c_L) 


( 11 . 21 ) 


= [x*,oo]^ + ^^/ dz^ z^[x^,z^]z{‘2.yF,j+iD^F,j){z^,z±)[z^,oo\-^ 


8ig 


as-^ 


as _ 

2 roo /•oo 


dz* / dz( z^x*, z^]zF,j{z^, z±)[z^, z'^]zF,^{z'^, z±)[zl,oo]z 


'Xf J Zt 

11.2.3 Gluon propagator in the shock-wave background 

The gluon propagator in a background gluon held (10.18) can be rewritten as 

1 


i{Al{x)A^{y)) = {x 


I \ab a:*>y* _ f ^ ia{x-y), 

P‘^ + 2igF + ie'^^f^'^ Jo 2a 


( 11 . 22 ) 


X (xx|e-*-^*Pexp{^ Hdz, e* - (A. + ^F) (z*)e-* - } 


'y* 


flU 


■ p\ , 

o^Pzy*\yj_yb 


Using the expansion (10.19) at y* = 0 we obtain with our accuracy 


1_=4^ ' 


F±^^{A.g^, + -F^,)e-^ 


P± 

as 


(11.23) 


~ 9y.p 


A. - —{2p>F,j - iD^F,j) 
as 


+ '^F^u + [2p>DjF^, - iDWjFf,, ) + ... 


Similarly to Eq. (10.21) we get 

1 


_|_ 2igF + ie 


2a \ 


^-ia{x-y). 


I \ab -m . 1 °° d:a ia' 

-\y)tiu = -i6{x^-y*) 1^ — + ie{y^-x^)l — 

X I d‘^z± ix±\e~'-^^*\z±){z±\g'^l{x^,y^]p±)FF^y*\y±){ll.24:) 


' —oo 


where 


Gyuix^,y^]p±) (11.25) 

f^* ( 2i ■ 

= 9yu[x*,y*]+ 9 dz* (- i^z^g^j,y{2p>[x^,z^]F,j{z^) - i[x^,z^]D^F,j{z^)) 

j y^ as 

2zz* 


H- 2^dip2v - dlp2y){[x^,z^]F,j{z^) H- -P^[x^, z^]DkF,j{z^)]Mz^,y^] 


as 


as 


g^2 pm* PZ:^: 2 

H- 3 / dz^ dz^ [ipf^uzi -P 2 /.P 21 /] K,^;*]F,j(z*)[z*,z(]F,^(z()[z(,y*] + 

Jy* ^y* 


Let us consider now quark terms coming from Fig. 7. From Eq. (10.24) it is clear that 
this contribution can be parameterized similarly to Eq. (10.31): 


(A“(x)7l^(y))Fig. 


’ da 


f-O 


da 


-6 »(x*-y*)/ —+0{y^-x^) — 

Jo 2a J_^ 2a J 


^-ia(x-y). / ^2^^ 


(11.26) 
(xx|e-*-"‘|zx) 


■Pj_ 

1 _i. -• 


2 2 2 . 

X (z±IQ'^t(x^,y*;p±)e^^^*ly±) + (y±le^^y*lz±)(z±IQf^(x^,y^;p±)e~^^^*lx±) 
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where and are given by expressions (10.32) - (10.37) with —)■ z* (and 

similarly —>• z'^). With our accuracy only the hrst terms in these expressions survive 

so and are given by Eq. (11.7) from previous Section. 

Adding gluon contribution (11.24) one obtains the hnal expression for gluon propagator 
in a shock-wave background: 




^-ia{x-y). 


(11.27) 


X / ci^zx 


2 2 

(x_L I I 2 ;_l ) ( 2 ;_l I (a:* ,y*-,p±) + (x*, y*; )] | y_L) 




■P± 
1 _ 


+ (y_L|e*-«^*|z_L)(2;_L|Q“^(x*,?/*;p_L)e *--^*|x_l) 


where is given by Eq. (11.25) and are given by Eq. (11.7). 

Similarly, for the complex conjugate amplitude one obtains (cf. Eq. (10.46)) 


da 


da 


{Al{x)A^{y)) = -e{y^-x^)l — + 0{x^ - y^) i — 


2a 


2a - 


^-ia{x-y). 


(11.28) 


X / d^z± 


.•d 




(xx|e-*-"* [a“^(x*,y*;px) + Q“td*,2/*;P±)]k±)d±|e'-"*|2/±) 


+ {y±\d^y*\z±){z±\ qI^ (x* , y*; )e"* I a;_L) 


where 


Gy.u{x*,y*-,p±) (11.29) 

= 9yA^*^y*] + 9 dz* [x*,2;*]j—^5'^i/(2F,jd*)d,y*]p^ + y*]) 

Jxt ^ Cl'S 

- 2^dlP2y - dip2y){F,j{zy,)[zy,,yy,] + 2i—DiF,j{zy,)[zy,,yy,]k‘-) [ 

OLS GLS 

8^2 ry* n, 2 ~ 

+ / dz^ / - -p2yP2iy] F,j{z^)[z^, zl]F,^{zi)[zi,y^] 


and are given by Eq. (11.10). Note that the transverse coordinates of all helds 

are effectively z±. 


11.2.4 Gluon emission vertex 

Similarly to Eq. (11.18) one obtains from Eq. (11.27) 


hm di(A“(A:)A^(y)) = y*, ^x;A:)(zx|e*d?/* |y^)e-s^- (11.30) 
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where 0“^(oo, y*, z_l; fe) is given by Eqs. (10.47)-(10.48). With our accuracy we get 


= g^,y[oo,y^] 


ab 


( r . . 4 ■ -1 \ 

+ ^ / dz^ ([oo,^*]^[-H- - Sip2f,)]F.j{z^, z±)[z^,y^]B 

I ni ''' CXS OlS / 

y* 

/ dz^ \ - ip2fiP2v[oo.,z^]zD^F,j{z^,z±)[z^,y^]z 
'y* '' 

f^* 4 I 

+ g dzl [2iag^„z'^ --p 2 fj.P 2 u\ [oo, z^]zF,j{z^, z±)[z^, zi]zF,^{z^, z±)[zi,y^]zj 


a^s'^ . 


as- 


^u: 


oo 


+ I dzl'ilj{z^,z±)[z^,oo]^F[oo,y^]^t''[y^,z'^]z'y^ f^i'y^'ilj{zl,z±) 


y* 'Jy* 

oo /•oo 


/•oo /•oo - 

- dz^ dz'^'ilj{z^,z±)'y^ ^i-f^[z^,y^]zt’'[y*,oo]X[oo,z'^]z'il^izl,z±)\ (11.31) 

J y* J z* ^ 

For the gluon emission in the complex conjugate amplitude one obtains (cf. Eq. (11.20) 

- lim k\{Al{x)Alik)) = 

fc2^.0 J 

(11.32) 

where 0“^(oo, y*, z_l;/ c) is given by Eqs. (10.49)-(10.50). With our accuracy 


C>‘iZix*,oo,z±;k) = fi(^j,[x*,(X)] 


ab 

z 


(11.33) 


/■“ / f 2iz^ _ . _ 

+ g dz^ f[x*,z*] 2 j— ^gfj,ui2y + iD^)F,jiz^, z±)[z^,oo]z 
Jxf ^ ^ Ck!'® 

- \idip 2 u - dlp 2 p){l + 2ik!-—bi)F,jiz^, z±_)[z^,oo\S\ 

OiS (y.s ^ ^ 

/ dz^ \^X^,Z^]z\^P2^lP2y&F,j{Z:^,Z2_)[z* 

Jxm ^ ^ 


ab 


4g 

a‘^s^ 


+ 


g[ dzi [-2iagf,uZ^ --p 2 fj.P 2 p] F,jiz^, z±)[z^, zijzPJ(zi, z±)[zi,oo]z\'] 
Jz» ® ^ ' 


ab 


aS‘ 


+ l^iy dz* I dzi'ipiz^,z±)[z^,x^]zt°'[x^,oo]zt’'[oo,zi]z'y-l;: ^ijlri’izl,z±) 


2 * 

OO /•OO 


dz^ / ciz(?^(z*,2;_L)7ll'l^i7M [2*,oo]2t^[oo,a:*]2t“[x*,4]2'0(2;(,2:_L 
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